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CHAPTER  1 
INTRODUCTION 

1.1     History  of  High  Temperature  Superconductivity 

The  year  of  1986  witnessed  one  of  the  most  revolutionary  discoveries  in  the 

history  of  science  -  high  temperature  superconductivity- by  Bednorz  and  Muller  [2] . 

The  critical  temperature  Tc,  below  which  the  electrons  can  pair  and  superconduct, 

was  unprecedentedly  as  high  as  35K  in  the  perovskite  structure.  Such  a  high 

transition  temperature  suggested  the  possibility  of  superconductivity  (SC)  even 

at  the  liquid  nitrogen  temperature  and  immediately  created  a  renaissance  in 

research  on  superconductivity.  The  enthusiasm  and  perseverance  invested  in  this 

area  over  the  past  two  decades  have  been  particularly  fruitful  and  have  led  to  the 

application  of  related  techniques  to  diverse  fields  of  industry.  Besides  the  previously 

inconceivable  increase  in  Tc  (  which  is  approximately  153K  now  !),  a  collection  of 

novel  physical  effects,  such  as  the  anisotropic  form  of  the  order  parameter,  peculiar 

normal  state  properties,  hidden  ordered  states  and  so  on,  emerged  successively, 

but  generally  in  unexpected  manners.  Their  appearances  kept  turning  over  our 

conventional  understanding  of  SC  within  the  BCS  framework  [3]  and  stimulated 

numerous  theoretical  concepts  aiming  at  the  explanation  on  the  underlying 

mechanism  of  HTS  [5,  6,  7,  8,  9],  However,  until  now,  a  basic  yet  convincing 

picture  of  HTS  has  not  been  found,  not  only  because  of  deficiencies  in  proposed 

theories,  but  also  because  of  the  lack  of  direct  experimental  tests  at  the  microscopic 

level  (Until  1999,  only  bulk  measurements  such  as  conductivity,  heat  capacity,  etc. 

were  available) .  Due  to  the  omnipresence  of  significant  disorder  in  all  HTS  samples, 

analysis  of  ah  such  experiments  required  a  model  of  disorder  treated  in  an  effective 

medium  approach.  This  obstacle  was  circumvented  in  1999,  by  the  notable  progress 

in  the  scanning  tunnelling  microscopy  technique.  The  differential  conductance  map 
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Figure  1—1:  The  crystal  structure  of  La2Cu04.  Red  sphere:  Cu2+;  blue  sphere: 
02~;  yellow  sphere:  La3+. 

displays  an  amazing  surface  structure  of  superconducting  samples  with  atomic 
resolution  [10].  The  improved  low  temperature  STM  boosted  our  observational 
capacity  to  an  unprecedentedly  fine  level  and  enabled  us  to  concentrate  on  the 
electronic  properties  of  high  temperature  superconductors  (HTS)  in  the  view 
of  local  probes.  This  has  raised  a  new  era  of  the  theoretical  study  of  HTS  and 
contributed  to  the  debate  over  the  nature  of  the  pseudogap  regime  where  electronic 
states  display  spontaneous  ordering  and  possibly  intrinsic  inhomogeneity. 
1.2     Background  About  High  Temperature  Superconductivity 
1.2.1     Crystal  Structure 

Among  the  various  types  of  HTS,  the  copper-oxide  compounds,  generally 
referred  to  as  cuprates,  are  of  particular  importance.  They  usually  possess  high 
critical  temperatures  and  complicated  crystal  structures.  But  what  really  makes 
them  so  interesting  is  the  illuminating  ubiquity  of  Cu02  planes.  These  layers  are 
separated  by  intervening  insulating  layers  (See  Fig.  1-1).  La2Cu04  is  a  typical 
candidate,  whose  structure  basically  mimics  a  ''sandwich,"  with  one  C11O2  plane 
between  two  LaO  planes.  Each  Cu  atom  is  surrounded  by  four  in-plane  and  two 
out-of-plane  oxygen  atoms  forming  a  perovskite  structure. 

The  electronic  structure  of  parent  compounds  of  HTS  is  quite  unusual.  Con- 
ventionally, transition  metal  oxides  have  'id  and  2p  orbitals  whose  bare  energies 
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are  well  separated,  and  the  strong  local  Coulomb  interaction  gives  rise  to  superex- 
change  mediating  long  range  antiferromagnetic  order  in  a  Mott-insulating  state 
with  gap  of  order  10  eV;  in  HTS  compounds,  Cu  atoms  lose  two  electrons  and 
become  divalent  ions,  leaving  a  hole  in  the  3dx2_y2  orbit  with  a  net  1/2  spin;  the  3d 
level  of  copper  ions  is  so  close  to  the  in-plane  2px„  level  of  oxygen  ions  that  they 
hybridize1    strongly  with  each  other  to  form  the  so-called  charge  transfer  insulators 
with  a  much  smaller  insulating  gap  (~  2ev). 
1,2,2     Phase  Diagram 

Superconducting  materials  are  formed  from  the  half-filled  antiferromagnetic 
Mott-type  ground  state  when  the  parent  compounds  are  doped,  leaving  mobile 
carriers  in  the  Cu02  plane.  The  evolution  of  the  electronic  and  magnetic  order 
upon  doping  establishes  a  complicated  phase  diagram  as  illustrated  in  Fig.  1 
2  [11].  The  intrinsic  antiferromagnetic  magnetic  order  is  quickly  destroyed  by 
adding  a  small  number  of  holes,  and  only  occupies  a  "sliver"  region  close  to  hall- 
filling,  whereas  antiferromagnetism  in  electron-doped  samples  is  more  robust. 
When  the  doping  increases,  a  spin-glass  behavior  may  be  present  before  the 
superconducting  phase  appears.  At  a  certain  minimal  hole  concentration,  we 
obtain  the  superconducting  state.  The  critical  temperature  increases  initially  with 
the  doping  concentration,  peaks  at  about  x  =  0.16  (optimally  doped),  and  then 
decreases  until  it  drops  to  zero.  Besides  the  high  critical  temperature  and  the 
peculiar  low  carrier  deusity  in  this  domain,  the  unusual  shape  of  the  energy  gap  in 
momentum  space  is  another  defining  property  in  HTS.  And  it  is  exactly  the  reason 


1  "Hybridize"  includes  many  physical  procedures  here:  the  crystal  field  split- 
ting lifts  the  degeneracy  of  3d  Cu  atomic  level;  Jahn- Teller  effect  further  lifts  the 
degeneracy  by  crystal  distortion;  the  resultant  3d  orbitals  then  hybridize  with  the 
in-plane  2p  oxygen  orbitals  to  form  Cu-O  bands  which  we  may  have  Mott  insulator 
or  charge  tmnsfer  insulator  depending  on  where  the  lower  anti-bonding  Hubbard 
band  lies. 
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Figure  1-2:  Schematic  plot  of  the  phase  diagram  of  hole-doped  high  temperature 
superconductors.  X  axis:  hole  dopant  concentration;  y  axis:  temperature. 

we  generally  refer  to  them  as  unconventional  superconductors:  the  order  parameter 
is  found  to  possess  a  dx2_y2  symmetry,  lower  than  that  of  the  Fermi  surface  itself. 

The  superconducting  phase  becomes  unstable  at  higher  temperatures  since 
thermal  fluctuations  break  the  ground  state  Cooper  pairs  and  a  "normal"  domain 
with  anomalous  attributes  emerges.  It  is  the  normal  phase  of  HTS  since  the 
temperature  T  >  Tc,  but  is  anomalous  in  the  sense  that  the  conventional  Fermi 
liquid  (FL)  theory  of  metals  fails  to  predict  its  one  particle  spectrum  and  transport 
properties  [12,  13,  14,  15].  Among  numerous  proposals,  the  so-called  "  Marginal 
Fermi  Liquid  theory"  [8]  (MFL)  was  aimed  to  address  the  normal  state  properties 
of  cuprates  by  reconciling  the  contradiction  between  the  well-defined  Fermi  surface 
and  the  non-Fermi  liquid  (but  universal)  behaviors  in  a  phenomenological  fashion. 
In  the  MFL  theory,  electrons  are  postulated  to  couple  with  bosonic  excitations 
whose  spectrum  is  only  a  scale-invariant  function  of  energy  and  temperature,  i.e., 
B(u,T)  ~  min(u>/T,  1).  As  the  consequence,  the  self-energy  for  the  coherent 
component  of  the  one  particle  spectrum,  extracted  from  this  proposal  displays  a 
peculiar  form  (See  Ch.  5  for  details)  with  logarithmically  vanishing  Z  —  (1  — 


5 
dReE((j)/dw)_1  at  zero  temperature  on  the  Fermi  surface.  The  Marginal  Fermi 
Liquid  theory  has  proved  to  be  reasonably  consistent  with  a  wealth  of  unusual 
experimental  observations,  such  as  the  linear  resistivity  [16],  in  several  decades  of 
temperature,  but  has  never  been  convincingly  derived  from  a  microscopic  theory. 
Especially,  how  it  yields  the  momentum-anisotropic  d-wave  pairing  symmetry  is 
the  intrinsic  deficiency  of  this  theory  (although  it  has  been  argued  that  this  could 
come  from  the  vertex  corrections) .  In  the  optimal-  and  overdoped  region,  while  the 
normal  state  property  is  fairly  depicted  by  MFL  theory,  a  conventional  BCS  theory 
is  believed  to  be  qualitatively  applicable  to  the  superconducting  states,  with  the 
order  parameter  taken  to  possess  d-wave  like  symmetry. 

There  is  one  well-established,  controversial  region  which  is  presently  at  the 
heart  of  the  debate  of  HTS,  which  lies  in  the  crossover  from  AF  order  to  SC  state, 
above  the  possible  messy  spin  glass  phase.  It  is  named  the  "pseudogap"  phase 
because  of  the  spectroscopic  identification  of  strongly  depleted  density  of  states 
(DOS)  near  the  Fermi  level.  The  pseudogap  temperature  T*,  below  which  this 
phenomenon  develops,  declines  linearly  with  increasing  doping  and  merges  with 
Tc  at  roughly  optimal  concentration.  This  phase  is  believed  by  some  physicists 
to  be  strongly  related  to  the  superconducting  state  in  that  it  evolves  smoothly 
into  the  latter  and  displays  a  similar  rf-wave  symmetry  in  the  order  parameter. 
So  far,  numerous  theories  such  as,  spin-fluctuation  [4] ,  resonant- valence  bond  [5] , 
precursor  scenario  [17],  and  competing  states  [18]  have  been  proposed  to  address 
the  underlying  mechanism  of  the  pseudogap  region  with  each  of  them  winning 
support  from  certain  experiments  [19].  Nevertheless,  there  is  yet  no  consensus 
about  the  origin  of  the  pseudogap  to  date. 

1.3     Current  Situation  in  High  Temperature  Superconductivity  Research 

The  research  on  HTS  in  the  past  twenty  years  attracted  an  enormous  amount 
of  theoretical  and  experimental  attention,  generating  interest  from  physicists, 
chemists,  material  scientists,  and  even  mathematicians  because  of  its  fascinating 
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Figure  1-3:  Current  situation  in  high  temperature  superconductivity. 

anomalous  behavior  observed  and  enticing  potential  applications.  However,  despite 
all  this  work,  no  consensus  on  the  origin  of  unconventional  superconductivity  has 
been  reached  and  there  is  a  huge  number  of  interpretations  on  HTS,  based  on  the 
diverse  observations  collected  from  different  groups,  with  different  techniques,  and 
on  different  physical  properties.  In  fact,  the  current  situation  in  the  study  of  high 
temperature  superconductivity  can  be  "figuratively"  depicted  by  Fig.  1-3  [20] , 
which  comes  from  a  famous  Indian  Buddhist  parable  about  the  blind  men  and 
the  elephant.  And  the  following  paragraph  from  this  famous  parable  reflects  how 
intensive  the  debate  among  them  over  the  shape  of  the  elephant  (  according  to 
everyone's  own  understanding)  is:   "Oh,  how  they  cling  and  wrangle,  some  who 
claim  For  preacher  and  monk  the  honored  name!  For,  quarrelling,  each  to  his  view 
they  cling.  Such  folk  see  only  one  side  of  a  thing." 

There  is  such  an  analogy  happening  in  the  high-Tc  community!  I  am  sure  that 
one  day  the  "Buddha"  (  a  correct  and  comprehensive  theory)  will  come  for  the  sake 
of  our  piety  and  opens  our  eyes  to  the  rest  of  the  "elephant"  (the  nature  of  high 
temperature  superconductivity),  but  I  would  like  to,  with  the  modest  expectation 
that  I  could  not  be  the  most  "blind"  disciple,  summarize  the  appearance  of  the 
monster  after  "caressing"  it  carefully: 
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1.3.1  What  We  Know  About  High  Temperature  Superconductivity 

1 .  Cuprates  have  copper  oxide  planes,  and  these  planes  are  crucial  to  high 
temperature  superconductivity; 

2.  The  electrons  pair  up  in  the  C11O2  plane  in  a  state  with  d-wave  symmetry; 
there  are  nodes  in  the  pairing  function  which  give  lower  energy  excitations  in 
quasiparticle  spectra; 

3.  The  carrier  density  in  HTS  is  considerably  lower  than  that  of  the  conven- 
tional superconductors;  the  coherence  length  is  as  small  as  £0  ~  3  —  4a; 

4.  It  is  the  proximity  to  a  underdoped  antiferromagnetic  Mott  insulator  that 
may  be  essential  to  understand  HTS,  and  probably  indicative  of  the  underly- 
ing magnetic  mechanism  of  HTS; 

5.  Cuprates  are  in  general  metals;  there  exists  generically  a  Fermi  surface  (at 
least  a  segment)  in  these  samples.  Normal  states  are  anomalous,  compared  to 
good  metals  which  can  be  well  described  by  Fermi  liquid  theory. 

1.3.2  What  Wo  Do  Not  Know  About  High  Temperature  Superconductivity 

1 .  What  is  the  mechanism  responsible  for  the  formation  of  Cooper  pairs? 

2.  What  is  the  nature  of  the  normal  state  and  the  pseudogap  regime? 

3.  What  do  the  disorders  in  HTS  tell  us?  What  is  the  source  of  the  local 
inhomogeneities  observed  in  the  tunnelling  experiments? 

The  collection  of  known  and  unknown  aspects  of  cuprates  paints  the  big  picture 
about  high  temperature  superconductivity  but  also  imposes  rather  constraints  on 
any  theory  regarding  the  fundamental  mechanism.  In  my  dissertation,  I  do  not 
intend  to  propose  any  novel  fundamental  mechanism  or  reconcile  the  discrepancies 
between  already  existing  theories;  instead,  I  will  concentrate  on  an  interesting  issue 
related  to  the  impurities  in  d-wave  superconductor,  which  has  been  the  focus  in  the 
recent  spectroscopic  experiments  on  cuprates  and  discussed  intensively  by  a  huge 
number  of  illuminating  works,  within  the  extended  version  of  conventional  BCS 
quasiparticle  theory  in  high  temperature  superconductivity. 


1.4      Disorder  in  Spectroscopic  Experiments 
1.4.1     Scanning  Tunnelling  Microscopy  (STM):  In  Touch  With  Atoms 

In  1981  Binnig  and  Rohrer  revolutionised  the  field  of  surface  science  by 
inventing  the  scanning  tunnelling  microscopy  [21].  When  a  fine  pointed  tip  of 
metallic  needle  is  brought  extremely  close  to  the  sample  surface  (in  general  a 
few  Angstroms)  and  a  voltage  is  applied  in  between,  the  wavelike  properties  of 
electrons  allow  quantum  tunnelling  to  happen.  The  probability  of  finding  the 
tunnelling  electrons  decays  exponentially  when  the  distance  between  the  tip  and 
the  sample  surface  increases.  The  topographical  image  of  the  surface  is  registered 
as  the  variation  in  the  current  when  the  tip  scans  over  the  surface.  First-order 
time-dependent  perturbation  theory  gives  out  the  tunnelling  current  as  function  of 
sample  bias  V  and  position  r: 

/(r,  V)  =  ^p  ^{/(M)[l-/(/J1,-el/)-y»[l-/(^-eK)]}|AV(r)|^(i?,-^+eV), 

(1.1) 

where  A/M1/(r)  is  the  tunnelling  matrix  element  and  /  the  Fermi  function.  The  5 
function  conserves  the  energy  as  it  does  in  the  usual  cases.  With  several  appropri- 
ate assumptions2  ,  we  can  write  /  as 

i-eV 

/(r.  V)  =  e-kz  I      dtp(r,  EF  +  e),  (1.2) 

Jo 

where  r  is  the  coordinate  of  points  sampled  and  z  is  the  vertical  distance  between 

the  tip  and  the  plane;  the  local  density  of  state  (LDOS)  is  defined  as 

p(rlE)  =  J2m*)\2S(E-E»)-  (1-3) 


2  1)  the  tunnelling  event  is  extremely  localized  (this  ensures  the  tunnelling  cur- 
rent is  proportional  to  the  local  density  of  states);  2)  it  is  spin- independent;  3)  for 
small  V,  we  can  using  WKB  approximation  such  that  A/(fc,  z)      =      Moe~*z,  where 
k2  =  4vnr(Aw  +  V)/h2  for  a  rectangular  barrier  of  height  A,„. 
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Assuming  that.  p(r,  E)  is  roughly  a  constant,  the  the  differential  tunnelling  con- 
ductance dljdV  is  proportional  to  the  LDOS  (this  can  be  obtained  by  moving  the 
tip  away  or  towards  to  the  surface  to  maintain  a  constant  current  and  then  apply 
an  alternating  voltage).  This  separation  enables  us  to  compare  the  theoretically 
calculate  LDOS  with  the  STM  measurement  qualitatively.  The  invention  of  STM 
technique  opened  an  new  era  of  surface  science  by  extending  our  "fingers"  to  atoms 
and  obtain  atomic  resolution  images  of  sample  surfaces  for  the  first  time.  It  was 
such  an  important  contribution  that  Binnig  and  Rohrer  were  awarded  with  the 
Nobel  Prize  of  physics  in  1986.  In  the  study  of  high  temperature  superconductivity, 
STM  also  contributed  in  a  unique  way,  especially  from  1999,  by  producing  stun- 
ning local  electronic  images  of  the  surface  of  BSCCO  and  providing  local  spectral 
information. 
1.4.2     Impurities  in  BiaS^CaCuaOs  :  A  Probe  of  High  Temperature  superconductivity 

Bi2Sr2CaCu208  is  a  typical  candidate  of  high  temperature  superconductors.  It 
is  well  known  that  the  Van  der  Waals  coupling  between  the  BiO  layers  in  BSSCO 
makes  it  easy  to  cleave  (Fig.  1-4)  and  obtain  the  image  of  its  surface  by  STM. 
However,  this  attribute  also  provides  space  where  dopant  oxygen  atoms  can  reside 
and  hence  result  in  interesting  experimental  observations.  Five  years  ago,  the  first 
high-resolution  STM  experiment  on  BSCCO  [10]  displayed  fascinating  atomic- 
scale  resonances  (as  shown  in  Fig.  1-5),  whose  spatial  and  spectroscopic  features 
agreed  crudely  with  the  theory  of  strong  quasiparticle  scattering  from  impurities 
in  a  d-wave  superconductor.  This  agreement  underscored  the  possibility  of  using 
impurities  of  different  electronic  properties  to  probe  the  superconducting  phase 
[22,  23].  Publications  about  STM  images  on  the  intentionally  impurity-doped 
BSCCO  soon  appeared,  supporting  the  simple  theoretical  impurity  scattering 
models  by  their  qualitative  agreement  with  predicted  resonance  energies  and 
fourfold  spatial  structure  [24,  25,  26].  Subsequent  STM  measurements  revealed 
further  the  existence  of  nanoscale  inhomogeneities  in  the  order  parameter  map 
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Figure  1  -4:  The  crystal  structure  of  BiaSrsCaCuaOg  . 

[27,  28,  29,  30]  which  are  currently  the  subject  of  intense  debate,  being  attributed 
either  to  interaction-driven  effects  such  as  stripe  formation  [31,  32]  or  to  the 
Priedel  oscillations  [28,  29]  of  weakly  interacting  quasiparticles  due  to  disorder. 
Those  results  were  then  Fourier  transformed  (FT-STS)  to  successfully  extract 
some  characteristic  wave-vectors  which  reflect  kinematics  of  quasiparticles  in  clean 
samples;  but  other  nondispersive  wave-vectors,  which  are  not  explicable  in  terms 
of  the  quasiparticle  scattering  scenario,  were  found  as  well  and  provoked  diverse 
proposals  concerning  the  ground  state  of  cuprates. 

Besides  STM  experiments,  angle-resolved  photoemission  (ARPES)  is  another 
leading  technique  in  the  study  of  HTS.  ARPES  provides  the  momentum  dependent 
single  particle  spectrum  and  information  regarding  the  self-energy  which  is  inherent 
from  many-body  effects.  In  the  ARPES  spectra  of  BSCCO,  the  width  of  the 
quasiparticle  peak  also  suggested  a  significant  elastic  component  in  the  self-energy 
due  to  the  scattering  from  spatially  extended  potentials  with  unknown  identity 
[33,  34]. 
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Figure  1-5:  STM  image  in  Bi2Sr2CaCu20s  showing  a)  atomic  scale  native  defect 
resonances;  b)  Zn  resonance.  Both  were  taken  at  zero  bias. 

Appreciable  differences  in  the  spectroscopic  and  transport  properties  between 
BSCCO  and  YBCO,  another  popular  HTS  material  (Tc  a  93A'),  were  reported 
extensively  [34,  35] ,  despite  the  fact  that  both  are  good  high  temperature  supercon- 
ductors with  well-identified  d-wave  order  parameters  as  well.  These  discrepancies 
were  not  successfully  reconciled  until  it  was  realized  that  unlike  YBCO,  BSCCO 
is  "intrinsically"  contaminated  by  disorder  through  the  way  it  is  prepared;  this 
was  gradually  established  by  the  accumulated  facts  from  STM  and  ARPES  exper- 
iments. The  charge  inhomogeneity  introduced  by  doping  is  poorly  screened  by  the 
states  in  the  C11O2  plane  and  induces  spatially  smooth  potentials  from  which  quasi- 
particles  scatter.  We  have  proposed  that  the  dopant  disorder  should  be  responsible 
for  the  observed  unusual  properties  of  BSCCO. 

Historically,  disorder  studies  have  played  essential  roles  in  the  comprehension 
of  superconductivity:  we  can  understand  the  mechanism  of  superconductivity  by 
understanding  how  it  could  be  destroyed  by  perturbations.  However,  until  recently, 
research  on  impurity  effects  was  limited  to  their  impact  on  bulk  properties  such  as 
transport  coefficients  or  heat  capacity,  in  which  a  finite  concentration  of  disorder 
is  generally  assumed  and  treated  in  a  statistical  manner.  It  is  the  improvement 
in  STM  image  resolution  that  enables  us  for  the  first  time  to  focus  on  the  local 
response  from  the  superconducting  state  to  one  single  impurity  individually,  from 
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which  we  can  attempt,  to  understand  inter-impurity  correlations  and  macroscopic 
disorder  phenomena  from  a  microscopic  point  of  view. 


CHAPTER  2 
SINGLE  IMPURITY  PROBLEM 

The  study  ol  the  impact  of  impurities  or  atomic-scale  disorder  on  the  physical 
properties  of  materials  has  been  of  great  theoretical  and  practical  significance.  In 
reality,  no  samples  are  pure.  As  Pauli  said:  "solid  state  physics  is  the  physics  of 
dirt."  Sometimes  impurities  are  destructive  to  desired  physics;  sometimes  they 
are  substituted  intentionally  to  bring  out  special  features  and  help  identifying 
the  underlying  physics  of  the  pure  system.  For  example,  impurity  induced  Friedel 
oscillations  can  aid  to  probe  the  kinematics  of  quasiparticles  in  metals  [36,  37] . 
In  this  chapter,  I  will  discuss  the  effect  of  a  single  impurity  in  an  otherwise 
clean  sample.  The  cases  in  normal  metals  and  conventional  superconductors  are 
reviewed,  followed  by  detailed  analysis  of  the  case  of  a  d-wave  superconductor.  I 
will  then  discuss  the  discrepancies  between  recent  STM  experiments  and  theoretical 
predictions. 

2.1     One  Impurity  in  Normal  State 
2.1.1     Bound  State  and  Resonant  State 

It  is  of  pedagogical  significance  to  review  the  problem  of  single  impurity  in 
a  metallic  state  before  we  begin  to  discuss  its  effect  in  a  superconductor.  The 
simplest  grand  canonical  hamiltonian  for  the  normal  state  problem  is 

H  =  £>(k)  -  ^{.A*  +  J2  Vocj^,,,  (2.1) 

k  k.k' 

where  cy,  ck  are  the  creation  and  annihilation  operators  of  particles  respectively. 
The  first  term  in  Eqn.  2.1  describes  the  single  particle  dispersion  of  the  non- 
interacting  ground  state  and  the  latter  denotes  the  on-site  energy  change  due  to 
the  perturbation  (The  single  impurity  is  assumed  to  be  located  at  R=0).  The 


13 


14 


Figure  2  1:  Diagrammatic  representation  for  the  scattering  process  by  single  im- 
purity. Coarse  arrow:  full  Green's  function;  solid  line:  pure  Green's  function;  cross: 
impurity;  dashline:  impurity  potential. 

single  particle  propagator  is  denned  as  G°(k,  k')  =  —  (Trcy{T)c'ki) ,  with  r  being  the 
imaginary  time.  In  the  case  of  free  fermions,  it  takes  the  simple  form: 

G°(k,«)  =  [w-e(k)r]  (2.2) 

The  single  impurity  problem  is  generally  solved  by  the  so-called  T-matrix  method 
which  provides  immediate  access  to  spectral  properties  [38].  The  perturbed  Green's 
function  can  be  expressed  in  terms  of  the  pure  propagator  G  and  the  T-matrix  as 

G  =  G°  +  G°TG°,  (2.3) 

where  the  T-matrix  is  defined  as 

T  =  V  +  VGT  (2.4) 

and  V  is  the  potential  matrix.  Fig.  2-1  shows  the  perturbative  expansion  of  the 
full  Green's  function  in  a  geometric  series  in  V. 

In  the  case  of  isotropic  point-like  scattering,  i.e.,  K(r)  =  Vo<S(r  -  r')  (Vo  is 
the  impurity  strength) ,  the  T-matrix  is  momentum-  independent  and  reduces  to  a 
simple  form: 

rM=i-voE%,,r  <2-5> 

The  poles  of  T-matrix  signify  the  formation  of  new  quantum  states,  since  they 
are  not  poles  of  Go,  i.e.,  the  eigenstates  of  the  pure  system.  There  are  two  kinds 
of  states  introduced  in  general,  as  shown  in  Fig.  2-2:  a  discrete,  real  bound 
or  anti-bound  state  located  outside  of  the  band  where  the  imaginary  part  of 
G"(u))  =  X]k(-'0(k''i')  vanishes,  i.e.,  where  the  density  of  states  (DOS)  is  zero  (the 
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Figure  2  2:  Schematic  illustration  ol  the  bound  state  in  normal  metal.  Leftisolution 
to  the  T-matrix;  Right:  the  bound  state  is  only  located  at  high  energies,  i.e.,  out- 
side of  the  band. 

right  panel),  and  an  overdamped  state  in  the  continuum,  which  has  a  finite  lifetime 
because  of  its  overlapping  with  the  background  (the  left,  panel) .  A  real  bound 
state  is  generically  localized  on  an  atomic  scale  if  the  impurity  potential  is  of  short 
range. 

2.1.2     Local  Density  of  States:  Friedel  Oscillations  and  Ripples  in  the  Fermi  Sea 
Apart  from  the  determination  of  bound  state  energies,  the  correction  to  the 
local  density  of  states  due  to  the  impurity  is  interesting  as  well.  The  introduction 
of  the  imperfection  does  break  the  stillness  of  the  degenerate  Fermi  sea,  inducing 
"ripples"  in  the  space  referred  to  as  the  Friedel  oscillations.  The  change  of  LDOS  in 
the  bound  state  is 

6N(t,w)  =  V02|Go(r,w)|2(5(w  -  Wo).  (2.6) 

In  a  two-dimensional  space,  the  real  space  Green's  function  is  evaluated  as 

G°(r,«)  =  -«rA'(0)//o1)[fc(w)r],fr(t<;)  =  kpjl  +  -.  (2.7) 

where  p,  is  the  chemical  potential  and  kf  is  the  Fermi  momentum. 

When  k(u)r  is  large,  the  Green's  function  has  the  following  asymptotic  form: 

2vr 


Cr(r,u) 


-iN(Q)[ 


k(u)r 


11/2    i[*Mr-»/4] 


(2.8) 


The  LDOS  falls  off  from  the  impurity  with  an  envelope  of  l/(kpr)  and  oscillates 
with  a  period  of  the  inverse  of  the  Fermi  vector. 
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Figure  2-3:  LDOS  modification  due  to  single  impurity  in  metals.  Right:  image 
on  the  surface  of  Ag  [36];  left:  theoretical  calculation  of  LDOS  around  a  single 
impurity  in  a  two-dimensional  metal. 

In  Fig.  2-3,  the  theoretical  prediction  on  the  distribution  of  LDOS  around  an 
repulsive  impurity  (right  panel)  at  zero  energy  in  a  two  dimensional  square  lattice 
is  shown.  The  four-fold  symmetry  of  LDOS  is  the  signature  of  the  underling  lattice 
symmetry.  In  the  left  panel,  the  Fourier  component  of  LDOS  (  obtained  from  the 
experiment  performed  by  Crommie  et  al. [36]  on  the  surface  of  Ag),  which  is  defined 
as 

p(q>u0<=5>,*p(r,w)  (2.9) 

r 

is  plotted:  the  radius  of  the  bright  feature  in  the  FT-STS  spectroscopy  is  2kp. 
2.2     Green's  function  for  BCS  Superconductors 
In  the  conventional  BCS  theory,  the  ground  state  of  superconducting  systems 
is  a  condensation  in  momentum  space  in  which  electrons  pair  up  into  "Cooper 
pairs."  The  mean  field  hamiltonian  for  the  BCS  superconducting  state  is  : 


HE 


Y,  £kcL<v + XKAi  ak  +  M. 


(2.10) 


with  At  =  ^(cidC-kj).  The  first  term  describes  the  motion  of  free  band  electrons 
and  the  second  term  describe  how  two  electrons  pair  up  to  form  the  superconduct- 
ing order. 
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It  is  convenient  to  define  the  column  vector  "Ek  and  its  conjugate1    as  the 
following: 

*k=        ""'     l*k  =  (<V-k,),  (2.11) 

and  work  within  this  framework. 

The  single  particle  Green's  function  is  then  a  2x2  matrix: 

G°(k,T)      =      -(TT<Pk(T)<S>l) 

_  I    {Ttcv(t)cI)       (TrCkT(r)c_kl)    \ 

\    (Trdki(T)Cv)       (TTclkl(T)C-ki)    j 

For  a  time-invariant  system,  we  can  farther  introduce  the  Fourier  transform: 


(2.12) 


G°(k,  ioj„)  =  J2  e'u'"rG0(k,  r),  (2.13) 

T 

where  lj„  is  the  Matsubara  frequency.  With  the  gap  Ajj  being  the  exact  self-energy 
due  to  the  phonos-mediated  electron-electron  interaction,  the  mean  field  single 
particle  Green's  function  for  the  clean  SC  state  (after  analytical  continuation  to  the 


real  axis)  is 


^«)=^^,  (2.14) 


where  the  quasiparticle  energy  defined  as  Et  =  ^/cjr  +  Ajr,  and  ro,ri,  73  are  the 
Pauli  matrices.  The  imaginary  part  of  the  real  space  Green's  function  (G(r,w)  = 


p(r^)  =  lYll>"^^  (2-15) 


1  This  spin  and  particle-hole  resolved  matrix  structure  is  conventionally  referred 
to  as  Nambu  representation. 


with  the  spin-resolved  LDOS, 

/>t(r,w) 


-7r-1ImG11(r,r,w  +  J0+) 
+7r"'lm  G22(r,  r,  -u  -  i0+) 


With  the  general  property  of  the  retarded  Green's  function: 

ImG(o)  +  i0+)  =  -haG(u  -  i0+) 

the  LDOS  may  also  be  written  as  : 
1 


P(r,w) 


7r_1Im  [Gu(r,  r,  u  +  i0+)  +  G22(r,  r,  -u  +  i0+)]. 
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(2.16a) 
(2.16b) 

(2.17) 

(2.18) 


2.3     Que  Impurity  in  an  s-wave  Superconductor 
2.3.1     Nonmagnetic  Impurity  Problem 

The  non-magnetic  disorder  problem  in  an  isotropic  s-wave  superconductor  was 
studied  by  Anderson  [39].  In  the  conventional  Nambu  notation,  the  hamiltonian 
with  the  presence  of  impurity  (at  R  =  0)  is  formulated  as 


H  =  Hb 


k.k'.a 


(2.19) 


where  Vo  is  the  strength  of  the  impurity.  It  was  pointed  out  that  in  dilute  concen- 
tration limit,  nonmagnetic  impurities  introduce  negligible  effects  on  bulk  properties 
(such  as  Tc)  of  isotropic  superconductors,  because  the  quasiparticle  energy  should 
be  conserved  in  elastic  scattering  and  the  momentum  transfer  of  quasiparticles 
in  scattering  process  cannot  disrupt  the  phase  coherence  established  if  the  pair- 
ing is  isotropic,  i.e.,  Ak  =  A0.  This  conclusion  goes  under  the  name  "Anderson 
Theorem."  As  a  trivial  test,  the  self-energy  of  quasiparticles  can  be  treated  with 
Abrikosov-Gorkov  formalism  [40]  (AG)  ,  yielding: 

1  I'kl 


T(k,u)- 


2tn  s/t'i  +  A* 


(2.20) 
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where  tn  is  the  normal  state  scattering  rate.  It  is  easy  to  see  that  the  quasiparticle 
elastic  scattering  rate  is  suppressed  from  its  normal  state  value  1/t/v  upon  the 
opening  of  a  gap,  and  it  is  zero  at  the  Fermi  level,  reflecting  Anderson's  theorem. 
2.3.2     Magnetic  Impurity  Problem 

When  a  magnetic  impurity  is  injected  into  the  superconducting  condensate, 
it  couples  with  the  spin  density  of  conduction  electrons  by  means  of  the  exchange 
interaction.  In  the  classical  limit,  2   the  impurity  term  in  the  momentum  space  may 
be  written  as: 

""•'  =  ^A?  E  J(k'  k')4»5W>  ■  Scic'/3,  (2.21) 

k,k' 

with  J(k,  k')  the  exchange  energy,  S  the  impurity  spin  and  a  the  spin  operator  of 
electrons.  Abrisokov  and  Gorkov  [40]  first  treated  the  magnetic  impurity  in  Born 
approximation,  followed  by  a  treatment  in  the  unitary-scattering  limit  by  Yu  [41] 
and  Shiba  [42] .  It  is  found  that  the  magnetic  moment  breaks  the  time-reversal  sym- 
metry upon  which  the  Cooper  pairs  are  formed  and  hence  lifts  the  spin  degeneracy 
of  quasiparticle  excitations.3    Consequently,  an  intragap  bound  state  appears  for 
spin  down  excitations  if  the  local  moment  is  spin  up  and  vice  visa.  Increasing  the 
concentration  of  magnetic  moments  leads  to  gapless  superconductivity  and  finally 
destroys  the  superconducting  state.  The  Tc  suppression  is  obtained  in  AG  form  [40] 


ln(srO  =  0(1/2)  -4>(\  +  pc),  (2.22) 


2  The  classical  limit  is  obtained  by  taking  J  — »  0  and  S  — >  oo,  but  keeping  JS  = 
constant.  In  that  sense,  the  spin  can  be  viewed  as  a  local  magnetic  field. 

3  This  can  be  verified  by  working  on  the  BCS  hamiltonian  in  the  Bogoliubov  for- 
malism. One  can  verify  that  the  magnetic  term  couples  in  the  Nambu  t0  channel 
and  changes  sign  for  spin  up  and  down  excitations.  A  potential  scatter  couples  in 
the  T3  channel,  introducing  nothing  but  level  shifts. 
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where  ij)  is  the  digamma  function'1  ,  Tc  and  Tm  are  the  actual  and  disorder- free 
transition  temperatures  respectively,  and  pc  =  (27rT(,Tc)_1  is  a  constant  determined 
by  the  self-energy  due  to  magnetic  interaction,  where  r^1  =  irTco/2y. 
2.4     One  Impurity  in  a  rf-wave  Superconductor 

Because  of  the  doping  procedures  through  which  superconducting  BSCCO 
samples  are  prepared,  atomic  scale  imperfections  are  inevitably  introduced. 
Recently,  the  differential  conductance  map  obtained  by  STM  experiments  on  the 
cuprates  BSCCO-2212  with  astonishing  resolution  have  verified  the  existence 
of  point-like  atomic  defects,  and  further  investigated  the  behavior  of  Zn  or  Ni 
atoms  in  the  d-wave  superconducting  bath  by  replacing  Cu  atoms  with  them 
on  purpose  [43,  44] .  I  will  mainly  review  the  works  on  single  impurity  in  d-wave 
superconductors  and  discuss  extensively  their  implications  for  STM  experiments. 
2.4.1     Nonmagnetic  Impurity 

Isotropic  Scattering  and  Pair-breaking  Effect.  It  is  speculated  that  when  a 
Zn  atom  (with  atomic  configuration  3d104s2)  replaces  the  on-site  Cu  atom  in  the 
copper  oxide  plane,  it  will  lose  2  electrons  and  the  resulting  Zn2+  should  behave 
like  a  classical  nonmagnetic  impurity  (if  we  neglect  its  quantum  nature).  This  is 
because  the  Zinc  cation  has  a  closed  3d10  shell  which  will  gives  a  zero  net  spin 
S  =  0.5 


4  The  digamma  function  is  defined  as 

0(«)    =    5jmr(z) 

oo 

=    -I  +  'V-t^ — r.*?4 -1.-2,-3 (2.23) 

where  7  is  the  Euler  constant. 

5  There  are  also  other  points  of  view  which  are  suggestive  of  the  magnetic  nature 
of  Zn2+  impurity  because  magnetic  moment  sensitive  experiments  [45]  revealed  that 
Zn  atoms  actually  induce  local  moments  and  hence  one  probably  has  to  also  take 
the  local  correlation  or  dynamic  effect  (Kondo  resonance)  into  consideration  [46]. 
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angle  on  FS    ty- 


Figure  2-4:  Schematic  illustration  of  the  pair  breaking  effect  due  to  potential  scat- 
ters in  a  d-wave  superconductor. 


However,  unlike  the  case  of  nonmagnetic  impurities  in  s-wave  superconductors, 

quasiparticle  scattering  due  to  Zn2+  ions  in  cuprates  violates  Anderson's  theorem 

without  breaking  time-reversal  symmetry.  The  physical  process  is  illustrated 

in  Fig.  2-4:  scattering  of  quasiparticles  mixes  the  initial  and  final  states  with 

different  momenta  k,  thus  mixing  order  parameters  with  different  phases.  When 

the  scattering  is  strong  enough,  it  mixes  the  states  of  complete  opposite  signs  of 

order  parameters  with  high  probability  and  suppresses  the  condensate.  However, 

1  emphasize  there  that  the  gap  structure  for  any  T  <  Tc  is  still  unrenormalized 

since  £)k  Ak  =  0  and  the  angular  structure  of  the  gap  function  remains  unchanged 

[47] .  To  illustrate  the  effect  of  impurities  in  superconductors  with  different  order 

parameter  symmetries,  I  hereby  tabulate  the  corresponding  Nambu  components  of 

self-energies  (in  Born  approximation,  i.e.,  S(w)  =  ntV^^kd"^^))  m  tne  taWe 

below: 

Table  2-1:  The  imaginary  part  of  self-energies  due  to  non-magnetic  impurities  in 
superconductors. 


So 

Ei 

IV      "      N 

0 

r,,     ° 

r„    A 

Va;-^ 

Vi'-i" 

where  T0  is  the  impurity  scattering  rate  in  the  normal  state,  with  Q  = 
u  -  Eo(w)  and  Ak  =  Ak  -  Sj.  The  symbol  (•)  denotes  angle  averaging  around  the 
Fermi  surface. 
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Gas  Model:  Circular  Fermi  Surface  and  Parabolic  Band.  The  single  impurity 
is  analytically  solvable  in  a  d-wave  superconductor  with  a  circular  Fermi  surface 
and  a  planar  continuum  form  of  the  gap  function,  i.e.,  At  =  Aocos(20).  With  the 
T-matrix  formalism  introduced  in  Sec  2.1.1  and  the  BCS  formalism  of  the  Green's 
function,  the  full  Green's  function  can  be  constructed  as 

G(k,k',w)  =  G0(k,Lo)6kk,  +  G°(k,u;)7>)G0(k\'-W.  (2.24) 

While  its  off-diagonal  component  vanishes  because  £]k  At  =  0,  the  T-matrix 
reduces  to  the  simple  form: 

T      =      T0T0  +  T3T3 

To    =    V02ga/(S+S-) 

T3    =    V02(e-g3)/(S+S„) 

So    =    ^E^V^  (2.25) 

k 

(2.26) 

where  l/(7riVoK>)  =  c  is  the  cotangent  of  the  s-wave  scattering  phase  shift  jjo, 
and  N0  is  the  density  of  states  at  the  Fermi  level.  Q  is  the  momentum  integrated 
Green's  function.  This  expression  has  resonances  when 


5±  =  1  -  V0(g3  T  So 


■■  0.  (2.27) 

0,  and  the  resonances 


In  the  special  case  of  a  particle  hole  symmetric  system  Q$ 
are  entirely  dominated  by  Qo(uj). 

In  the  unitary  limit,  i.e.,  %  -»  7r/2  (this  is  equal  to  c  <  1),  the  resonance 
energy  f!„  and  scattering  rate  T  are 

±7TCA0 


r 


21og(8/«) 

tAAq 
41og2(8/7rc)' 


(2.28a) 
(2.28b) 
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These  results  were  first  obtained  by  Balatsky  et  al.  [24].  Clearly,  the  two  symmet- 
ric resonances  are  the  signature  of  electron-hole  duality  nature  of  BCS  quasipar- 
ticles;  furthermore,  there  is  a  finite  damping  rate  for  those  resonances  since  they 
emerge  from  a  continuum  background,  namely  p(u)  ^  0.  However,  a  salient  feature 
of  those  resonances  should  be  pointed  out:  when  c  — >  0,  the  resonances  will  be 
tuned  towards  to  the  Fermi  level  and  the  damping  rate  diminishes  simultaneously, 
i.e,  the  resonant  states  become  "marginally"  bound.  These  rnidgap  states  are 
generally  referred  as  virtual  bound  states  because  of  their  asymptotic  behavior 
described  above,  which  can  be  attributed  to  the  linearly-vanishing  DOS  (p(u>)  ~  u) 
at  low  energies. 

The  spatial  distribution  of  the  impurity  induced  states  conveys  additional 
information  about  the  nature  of  HTS  and  can  be  measured  directly  by  STM 
experiments.  The  theory  anticipates  that  the  LDOS  spreads  in  an  anisotropic 
manner  [24,  26] :  it  decays  as  the  inverse  second  power  of  the  distance  from 
impurity  along  the  nodes  of  the  gap  and  exponentially  in  the  vicinity  of  the 
extrema  of  the  gap,  and  the  decay  length  of  the  latter  is  characterized  by  £0  = 
vF/fiA0.  The  spatial  patterns  of  these  "native"  defects  [10]  agree  with  the  naive 
theory  to  the  extent  that  the  observance  of  a  bright  spot  of  the  diameter  of  2- 
3nm6  ,  with  p(r)  decaying  s  1/r2  for  r  >  £0,  i-e.,  the  Friedel  oscillations  at  the 
periphery  of  the  central  bright  image,  as  illustrated  in  Fig.  2-5. 


The  coherence  length  extracted  from  the  size  of  the  atomic  resonance  is  around 
15  A.  This  again  coincided  with  the  £o  known  for  Bi2Sr2CaCu208  ,  being  a  convinc- 
ing evidence  for  the  quasiparticle  scattering  scenario.  However,  the  anisotropic  con- 
figuration of  LDOS  distribution  within  a  distance  of  {0  from  the  impurity,  namely, 
the  structure  inside  the  bright  spots,  was  not  discernible  due  to  the  poor  spatial 
resolution.  Later  improvement  in  STM  yielded  differential  conductance  map  with 
subatomic  details,  and  the  fine  structure  of  LDOS  around  impurities  were  obtained 
eventually,  which  invoked  another  quarrels  upon  the  tunnelling  mechanism  and  the 
nature  of  impurity  itself. 
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Distance  from  scattering  center  |nm) 

Figure  2-5:  The  fall-off  of  LDOS  around  a  native  defect. 

To  study  the  spatial  modulation  of  LDOS,  we  need  the  explicit  form  of  the 
Green's  function.  These  have  been  worked  out  by  Joynt  and  Balatsky  in  some 
special  cases  [26,  48].  I  also  evaluated  the  Greens  function  when  for  the  situation 
ui/Aq  <  1/kpr  <  l/fef^o-  The  results  are  tabulated  as  follows  for  future  adoption: 


G°(R,u)  (2.29) 

n°t£z$%$i*  R II  ("0) 

""7^  [(«  A^7r°  +  r,  +  T3)  cos  kFR+    R  ||  (100) 
(''Afr7"0  +  ri  ~  T3>  sin^F-Rj 
Lattice  Model,  Strong  scatters  and  Unitary  Limit.  The  copper-oxide  plane 
of  cuprates  materials  consists  one  Cu  atom  and  two  oxygen  atoms  in  one  unit 
cell.  A  square  lattice  with  copper  atoms  on  the  commensurate  sites  is  the  simplest 
reasonable  model  for  practical  computations.  While  the  real-space  components 
of  the  Green's  functions  are  analytically  unobtainable  in  the  SC  state,  I  hereby 
present  the  numerics  for  both  a  simple  tight  binding  model  ek  =  -2t(coskT  + 
cos  ky)  -  n  and  a  realistic  6-parameter  tight-binding  model  proposed  by  Norman  et 
al.  [49],  both  having  the  corresponding  d-wave  order  parameter  Ak  =  A0(cosfcI  - 
cosfcv).  Note  the  maximum  value  of  the  order  parameter  in  the  half-filled  lattice 
system  with  the  current  convention  is  2A0. 
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Figure  2-6:  On-site  and  NN  site  spectra  for  tight  binding  band  with  one  impurity. 

Fig.  2-6  shows  the  results  for  the  nearest  neighbor  tight  binding  model.  The 
impurity  potential  is  taken  to  be  V0  =  5.3(i,  which  generates  two  resonances  at 
^     =  ±0.013t]  ~  ±1.5men  using  the  Norman  band,  chosen  to  reproduce  the 
observed  resonance  energy  of  Zn.  The  LDOS  on  the  impurity  site  and  its  neighbor 
site  are  plotted.  Note  that,  the  weight  is  almost  excluded  from  the  impurity  site  due 
to  the  strong  potential  7   and  the  peak  at  Q  =  -0.013d!  captures  the  essence  of 
resonant  physics.  Because  the  off-diagonal  components  of  the  integrated  Green's 
function  vanish,  i.e.,  Gi2(r  =  0,u)  =  0,  only  the  11  element  of  the  T-matrix 
contributes  to  the  LDOS  and  this  yields  only  one  peak  in  the  on-site  spectrum. 
The  situation  is  completely  inverted  on  its  nearest  neighbor  sites,  where  one 
observes  a  pronounced  peak  at  positive  sample  bias.  The  coherence  peak  heights 
are  also  strongly  suppressed,  implying  that  superconductivity  is  locally  suppressed 
as  well.  There  are  two  distinct  ranges  for  the  falloff  of  LDOS:  an  exponential 


7  The  integrated  spectrum  weight  on  arbitrary  site  is  still  conserved 
/  dwp(r,uj)     =     1.  The  sum  rule  is  verified  by  the  excessive  weight  trapped  at  high 
energy  regime,  i.e.,  the  real  bound  or  antibound  states.  However.we  focus  on  the 
low  energy  bound  states  within  the  gap. 
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envelope  along  the  antinodal  direction  and  a  power  law  decay  along  the  gap 

nodes.  Additionally,  the  alignment  of  LDOS  around  the  scatters  shows  a  four-fold 

symmetry. 

2.4.2     Magnetic  Impurity 

Unlike  the  nonmagnetic  impurities,  even  for  a  s-wave  superconductor,  mag- 
netic impurities  are  pairbreaking  since  they  violate  time-reversal  symmetry,  lift  the 
Kramers  degeneracy  of  the  quasiparticle  states  and  hence  introduce  intragap  bound 
states.8    While  the  theoretical  expectation  for  a  magnetic  impurity  in  a  s-wave 
superconductor  (for  example,  Nb)  has  been  experimentally  verified  by  Yazdani  et 
al.  [50],  Hudson  et  al.  [44]  performed  the  first  atomically  resolved  STM  experiment 
on  Ni  doped  Bi2Sr2CaCu208  (rf-wave  superconductor).  As  shown  in  Fig.  2-7,  two 
distinguishable  peaks  were  reported,  at  fii  =  9.2  ±  1.1  meV  and  fi2  =  18.6  ±  0.7 
meV  separately  in  both  the  particle  and  hole  channels.  The  spectral  weight  of  the 
impurity  state  decayed  in  an  oscillating  manner:  these  peaks  were  particle-like  (i.e., 
on  positive  sample  bias)  on  the  impurity  site,  then  became  hole-like  on  the  next 
nearest  sites  and  again  particle-like  at  the  2-NN  sites.  Therefore,  the  LDOS  around 
the  impurity  at  positive  and  negative  biases  were  rotated  by  45°  with  respect  to 
each  other.  The  spatially  complementary  feature  led  to  overall  nearly  particle-hole 
symmetric  site-averaged  spectra;  moreover,  the  coherence  peaks  remained  robust. 
A  conclusion  which  was  apparently  in  conflict  with  our  naive  intuition  was  then 
drawn:  superconductivity  was  not  locally  disrupted  substantially  by  the  Ni  impu- 
rity. Detailed  analysis  even  revealed  that  the  potential  scattering  is  the  dominant 
agent  in  Ni  doped  cuprates  since  the  associated  energy  is  of  roughly  5-10  times 
bigger  than  those  of  other  channels. 


8  We  neglect  the  Kondo  physics  here.  The  Kondo  screening  in  HTS  supercon- 
ducting state  is  novel  in  sense  that  it  takes  place  on  the  background  with  depleted 
DOS  rathern  than  on  the  background  with  a  roughly  constant  DOS  around  the 
Fermi  level. 
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Figure  2  7:  Conductance  spectra  above  the  Ni  atom  and  at  several  nearby  loca- 
tions. Solid  circle:  on-site;  open  circle:  NN-site;  square:  2NN-site;  triangle:  30A 
away  from  Ni;  solid  line:  average  of  the  first  four  spectra. 
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Figure  2-8:  Left:  the  schematic  plot  of  Cu02  plane;  middle:  experimental  LDOS 
around  Zn  impurity;  right:  theoretically  predicted  LDOS  for  V0    =    5.3t  (rotated  by 
45"  with  respect  to  a)  and  b)). 

2,4.3     Discussion:  Characteristics  of  Impurities 

Classic  impurities  are  generally  categorized  into  non-magnetic  and  magnetic 
ones,  with  the  Zn  and  Ni  as  the  representatives.  Although  experiments  exploiting 
them  as  probes  for  the  superconducting  state  have  achieved  crude  agreement 
with  the  quasiparticle  scattering  theory,  characterization  of  their  own  nature 
is  far  from  complete  yet.  While  the  measurements  on  the  bulk  properties  such 
as  resistivity  [51,  52],  microwave  conductivity  [53]  and  Tc  suppression  display 
qualitatively  similar  dependences  of  Zn  and  Ni,  the  magnetically  sensitive  probes 
such  as  NMR  (nuclear  magnetic  resonance)  [45],  INS  (inelastic  neutron  scattering) 
exhibit  appreciable  differences  between  Zn-  and  Ni-doped  samples.  On  one  hand, 
the  standard  Abrikosov-Gorkov  formalism  predicts  that  Tc  is  only  suppressed  about 
20%  faster  by  Zn  than  Ni  if  the  magnetic  nature  of  Ni  impurity  is  neglected;  on 
the  other  hand,  the  superfluid  density  is  known  to  be  strongly  depleted  by  Zn 
rather  than  Ni.  The  image  of  the  immediate  surroundings  of  Zn  and  Ni  provided 
by  the  STM  experiments  now  provides  a  coherent  microscopic  picture:  while  Zn 
is  fatal  to  local  superconductivity,  Ni  has  a  much  weaker  impact,  consistent  with 
earlier  transport  measurements  suggesting  Zn  is  a  unitary  scatterer  and  Ni  is  an 
intermediate  strength  scatterer. 

Besides  those  differences  between  Zn  and  Ni  imparities,  the  sign  of  Zn  poten- 
tial itself  is  a  matter  of  controversy.  It  was  first  claimed  that,  according  to  Hudson 
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et  al,  Zn  atoms  are  attractive  potentials  for  electrons  since  "the  on-site  (impurity 
site)  spectrum  only  has  one  peak  at  the  negative  sample  bias  around  fi0  =  -1.5 
meV  ±0.5  nieV  "  [10].  However,  theoretical  calculations  predict  that  an  attrac- 
tive potential  should  break  the  pairs  and  bind  the  electrons,  yielding  prominent 
tunnelling  DOS  peaks  at.  positive  sample  bias  [24,  54].  Beyond  the  position  of 
the  resonance  peaks,  there  are  other  discrepancies  between  the  STM  differential 
conductance  map  and  the  predicted  LDOS.  In  Fig.  2-8,  the  schematic  plot  of  the 
structure  of  CuOj  plane  is  given,  where  the  orange  spots  and  green  spots  stand 
for  copper  atoms  and  oxygen  atoms  respectively.  The  middle  panel  shows  the 
experimental  results  around  a  Zn  impurity  with  the  same  lattice  orientation  as  the 
left  panel  and  the  rightmost  panel  displays  the  theoretical  calculation  of  LDOS 
for  a  repulsive  potential  (Va  =  5.35t),  but  rotated  by  45°  degrees  with  respect 
the  first  two  panels.  Strong  potentials  scatterers  are  expected  to  exclude  weight 
significantly  from  the  impurity  site  and  produce  an  on-site  LDOS  minima  with  a 
relatively  bright  spots  on  the  next  nearest  neighbor  sites,  as  illustrated  by  the  right 
panel  in  Fig.  2-8,  whereas  the  experiments  observe  completely  reversed  patterns 
with  a  local  maxima  on  the  impurity  site  and  darker  spots  on  the  NN  sites  (the 
middle  channel) .  Comparison  between  the  middle  and  right  panels  elucidates  these 
discrepancies  clearly. 

Is  Zn  a  repulsive  or  attractive  potential  in  BSCCO-2212?  A  hand-waving 
argument  can  be  given  by  considering  the  atomic  configurations  of  Zn  and  Cu 
atoms.  After  losing  two  electrons,  the  energy  of  the  closed  'id  shell  of  Zn  should 
be  a  priori  far  below  the  Fermi  level,  as  derived  from  the  tabulated  ionization 
energies  (Table  2-2)  of  Copper  and  Zinc  atoms,  since  the  energy  to  create  a  Cu++ 
is  less  than  that  required  to  create  Zn+4 .  It  is  then  speculated  that  Zn  atoms 
should  be  attractive  potentials  relative  to  background  Cv?+  ions  with  strength  of 
approximately  -2ev. 
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Table  2-2:  The  ionization  energy  of  Copper  and  Zinc  atoms. 

Cu  Zn 

Atomic  configuration  3dw4sl  3rflu4s2 

first  ionization  euergy(eV)         7.720  9.934 

second  ionization  energy(eV)     20.292  17.904 

third  ionization  energy(eV)        36.83  39.722 

However,  calculations  based  on  the  criterion  for  unitarity  in  the  realistic  band 
structure  of  BSCCO  [49]  yield  an  opposite  result.  It  has  been  well  established  by 
experiments  that  Zn  atoms  are  indeed  hard-core,  unitary  scatters,  for  example,  the 
phase  shift  T)0  extracted  from  the  STM  data  of  Zn  impurity  [44]  is  0.487T,  close  to 
7r/2,  i.e.,  unitary  limit.  One  can  also  perform  computations  on  the  phase  shift  with 
the  conventional  T-matrix  fashion 

,    ,  ,  Im  detf 

'/M  =  tan"1  ..  (2.30) 

Re  detT 

It  is  found  that  only  a  repulsive  potential  can  possibly  approach  to  this  limit.9 

In  Fig.  2-9,  I  show  the  resonance  energy  as  a  function  of  the  impurity  potential, 

as  the  resonance  energy  is  defined  to  be  the  position  of  the  local  maxima  in  the 

on-site  (the  impinity  site)  spectrum.  It  is  clear  that  to  drive  the  resonant  frequency 

towards  to  the  Fermi  level,  a  repulsive  potential  with  V0  ss  20|<i|  is  needed.  Finally, 

a  recent  ab  initio  calculation  exploiting  density  functional  theory  on  the  band 

structures  and  effective  potentials  of  different  defect-doped  BSCCO  samples  [55] 

claims  that  Zn  atoms  are  actually  short  range  repulsive  scattering  centers. 

Recently,  two  similar  arguments  [56,  57]  were  addressed  to  reconcile  the 

discrepancy  between  the  result  of  attractive  potentials  and  the  spectra  from 

experiments.  It  was  contemplated  that  the  spectral  weight  on  the  nearest  neighbor 

Cu  sites  in  the  CuOj  planes  appears  above  the  surface  Bi  atoms  directly  over  the 


9  The  conclusion  is  rigorously  robust  regardless  of  the  detailed  band  structure 
as  long  as  the  particle-hole  asymmetry  is  introduced  and  a  hole-doped  sample  is 
discussed. 
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Figure  2-9:  The  resonance  frequency  as  a  function  of  scattering  potential  V0;  inset: 
the  resonance  frequency  as  a  function  of  1/V'0. 

impurity  sites  since  STM  tips  are  indeed  probing  the  wave  functions  or  LDOS  of 
the  nearest  sites  rather  than  the  impurity  site  itself,  the  so  called  "filter  effect." 
The  Zn  impurity  is  then  identified  as  an  attractive  scatter,  since  it  is  the  "filter 
effect"  that  allows  us  to  observe  a  pronounced  on-site  peak  at  negative  energy 
(an  attractive  potential  is  expected  to  yield  a  peak  on  the  positive  sample  bias) 
However,  it  doesn't  explain  why  the  calculations  based  on  the  realistic  band 
structure  aiming  to  the  unitary  limit  necessitate  a  repulsive  potential. 


CHAPTER  3 
QUANTUM  INTERFERENCE  BETWEEN  IMPURITIES 

Images  of  the  LDOS  around  impurities  have  confirmed  the  existence  of 

resonant  quasiparticle  states  near  strong  scatters  as  Zn,  but  have  given  rise  to 

new  questions  regarding  the  microscopic  model  for  impurities  as  well.  While  the 

debate  on  single  impurity  is  unsettled  yet,  HTS  materials  usually  contain  finite 

concentrations  of  impurities  and  the  impurity  wavefunctions  do  interfere  with  each 

other.  The  most  remarkable  consequence  of  the  quantum  interference  between 

many  impurities  is  the  formation  of  a  impurity  band  at  low  energies  (the  nodal 

quasiparticle  states),  which  is  believed  to  dominate  the  transport  properties  such 

as  microwave  conductivity.  The  many-impurity  problem  has  been  treated  with 

numerous  methods,  including  the  self-consistent  T-matrix  approximation  (SCTMA) 

[48,  58,  59,  60],  exact  diagoualization  of  Bogoliubov  De-Gennes  equations  (BDG) 

[61,  62,  63],  nonperturbative  techniques  [1,  64,  65]  and  perturbative  treatment 

including  weak  localization  effects  [66].  None  of  these  approaches  offers  insight  into 

the  mechanism  of  interference  itself  and  we  are  still  unable  to  answer  the  question 

of  how  the  macroscopic  disorder-averaged  properties  of  a  d-wave  superconductor 

actually  arise  from  the  local  properties  of  individual  impurity  states.  The  lack  of 

the  connection  between  the  understandings  of  physics  at  these  two  different  scales 

initiated  the  studies  on  the  two-impurity  problem  [54,  67,  68],  which  displays  the 

quantum  interference  effect  at  the  simplest  level.  In  this  chapter,  I  will  report  my 

study  on  the  two-impurity  problem  with  a  nongeneric  band.  I  will  then  generalize 

the  conclusion  from  this  special  example  to  generic  cases  and  discuss  how  we  can 

make  the  bridge  connecting  the  single  impurity  result  to  the  STM  experiments  on 

the  fully-disordered  systems. 
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3.1     Two  Impurities  in  a  rf-wave  Superconductor 
3.1.1      r-matrix  and  Resonance  Frequency 

With  the  formalism  introduced  in  the  previous  chapter,  the  T-matrix  for 
two  isolated  point-like  impurities  can  be  obtained  by  iterating  the  procedure  for 
single-impurity  T-matrix,  in  a  4  x  4  basis  of  spin  and  impurity  site  labels, 

rf,    ,   ,       I  ff'  fTcG°(K)f„,  \ 

Ttm{w)=  \  (3.1) 

\^  fTmG°(-R)Te  Jfm         j 

where  R  =  Rt  -  R„  and  where  Te,  Tm  are  the  single  impurity  T-matrices  associated 
with  the  two  impurities.  For  identical  impurities,  ft  =  fm  =  f(u),  the  single 
impurity  T-matrix  defined  previously.  The  quantity  /  is  defined  as: 

/(«)  =  [1  -  G°(-R,u;)fiMG0(R,W)f,„(w)]-1,  (3.2) 

where  G°(R,w)  =  £kexp[ik  '  R]GkM  is  .i|lst  tne  Fourier  transformation  of 
G^(u),  the  unperturbed  Nambu  Green's  function.  For  systems  with  inversion 
symmetry  G°(R,ui)  =  G0(-R,w).  Note  that  in  Eqn.  3.1,  the  physical  processes  are 
clearly  identifiable  as  multiple  scatterings  from  each  impurity  I  and  m  individually, 
plus  interference  terms  where  electrons  scatter  many  times  between  I  and  m.  In 
k-space,  we  can  write  the  T-matrix  in  the  more  usual  2x2  notation  as 

TWH  =  [e'k-R'r0    e'kR'"r0]T>m       '  '  (3.3) 

e-.k<.R,„ro 

where  to  is  the  Pauli  matrix. 

Provided  the  resonance  energies  are  distinct,  peaks  in  the  total  density  of 
states  correspond  to  minima  of  the  T-matrix  denominator: 

V  =  det[l  -  G°(-R,uj)f  (u>)G°{R,u)f  («)].  (3.4) 


Explicitly,  V  =  Z>i2V(S+Sl )  with 


where 


£>i     =    V<V;  +  V02G'iCR,uj) 

V±    =    [l-KoGfs(0,w)±VoGo(0,w)] 

+(-1)"Vo[=fG0(R,o;)  +  G3(R,u)]. 
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(3.5) 


(3.6) 


The  factors  D] ,  T>2  determine  the  four  2-impurity  resonant  energies. 

Here  Ga(R,u)  is  the  r„  component  of  the  integrated  bare  Green's  function 


GQ(R,w)  =  -Tr(raG0(R,w)). 


(3.7) 


In  certain  special  configurations,  e.g.  if  the  two  impurities  are  located  at 
45°  with  respect  to  one  other,  it  is  easy  to  check  that  the  off-diagonal  Green's 
function  Gi(R,w)  =  0  V  R.  In  this  case  the  entire  resonant  denominator  factories 
V  =  Z>i+Z>i_Z>2+©2--  The  T-matrix  then  takes  the  simple  diagonal  form 


Pl- 

-  + 

T+ 

T_ 

V2- 

+ 

&H. 

fk,k-M    =   2V„cos(k-5:)cos(k'.5;) 

+    2K„sin(k.5:)sin(k'.5:) 

(3.8) 

where  r±  =  (r3  ±  r0)/2. 

When  two  identical  impurities  with  resonance  energies  fio-^o  are  brought 
together,  the  bound  state  wavefunctions  interfere  with  one  another,  in  general  split- 
ting and  shifting  each  resonance,  leading  to  four  resonant  frequencies  Qj",  Qj",  fij 
and  Q.| ,  where  the  subscript,  indicates  which  factor  in  Eqn.  3.6  is  resonant.  If 
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splittings  are  not  too  large,  the  electron  and  hole  resonances  are  related  in  a  sim- 
ilar way  as  in  the  1-impnrity  problem,  fif  ~  -fij  arid  Q2   ~  -Qj .  Again  the 
weight  of  each  resonance  may  be  quite  different  or  even  zero  on  any  given  site.  A 
large  splitting  may  be  taken  as  evidence  for  strong  hybridization  of  quasiparticle 
wavefunctions.  If  we  take  the  interimpurity  distance  R  as  a  parameter  and  keep 
impurity  potentials  and  other  parameters  fixed,  there  are  two  obvious  limits  where 
this  splitting  vanishes.  In  the  case  of  separation  R  =  0,  the  two  impurities  combine 
(mathematically)  to  create  a  single  impurity  of  strength  2V0,  so  both  Qf2  approach 
the  Qg  (2Vo)  appropriate  for  the  double  strength  potential.  In  the  case  of  infinite 
separation  R  — i  oo,  we  must  find  Q^2  approaching  the  fio(Vb)  appropriate  for 
isolated  single  impurities.  To  illustrate  this  argument  explicitly,  I  will  present  both 
the  analytical  calculation  with  a  gas  model  and  numerical  results  on  a  tight  binding 
lattice  model. 
3.1.2     Resonance  Frequencies 

Gas  Model.  Eqn.  3.4  is  a  general  result  for  two  ^-function  potentials  embedded 
in  a  host  described  by  an  arbitrary  G°.  We  would  like  to  derive  analytical  results 
for  the  resonance  energies  obtained  therefrom  to  get  some  sense  of  the  appropriate 
length  scales  and  symmetries  in  the  problem.  At  large  distances,  the  resonance 
energies  must  approach  the  single  impurity  values,  so  the  splittings  can  be  cal- 
culated perturbatively.  To  do  so  one  must  first  obtain  analytical  expressions  for 
the  large-distance  behavior  of  the  unperturbed  Green's  functions.  This  is  difficult 
for  the  superconducting  lattice  tight-binding  model  on  which  most  of  this  work  is 
based,  but  much  insight  can  be  gained  by  studying  the  equivalent  gas  model,  with 
spectrum  ek  =  k2/2m.  In  this  case  expressions  have  been  obtained  by  Joynt  and 
Balatsky  et  al.  [24,  48]  for  the  d-wave  integrated  Green's  functions  Ga(R,ui  =  0) 
at  large  distances,  both  for  R.  making  an  angle  45°  or  0°  with  the  x  axis.  We  have 
evaluated  the  real  space  Greens  function  in  the  last  chapter  for  the  range  of  fre- 
quencies w/Ao  <  1/fcf-r  <C  1/ftffo-  With  these  expressions,  the  resonance  energies 
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may  now  be  found  by  inserting  these  expressions  for  frequencies  u  =  Q*  +  S  into 
(3.4)  and  solving  for  the  shifts  S.  We  find  fi^  ~  fij  ±  S,  with 

(3.9) 

'    rj[      w^r^+t^       RD("0) 

A'""J'     1  vwZ^,)  «'s(fc^  +  */4)  R II  (ioo) 

These  expressions  are  valid  for  S/Uq  <g  1. 

Clearly  the  decay  of  the  splitting  ~  exp(-r/£o)/v,**»"  is  much  more  rapid  for 
distances  larger  than  the  coherence  length  along  the  antinode  (100)  than  for  along 
the  nodes,  where  it  falls  as  ~  1/r.  Tliis  anisotropic  form  of  the  splitting  is  the 
signature  of  the  marginal  bound  (nearly  localized)  single  impurity  state  in  a  rf-wave 
superconductor,  which,  by  allowing  propagating  low-energy  nodal  quasiparticles, 
has  extended  tails  along  the  gap  nodes.  It  is  the  lack  of  a  scale  along  this  particular 
direction  in  the  impurity  wavefunctions  that  allows  the  strong  overlapping  even 
over  long-distances.  We  will  see  later  that  how  this  attribute  and  its  consequence 
questions  the  STM  analysis  of  "isolated"  impurity  patterns. 

Lattice  Model.  In  the  following  sections,  I  will  briefly  sum  up  the  numerical 
results  on  two  types  of  lattice  band  structure:  the  nearest  neighbor  hopping 
model  (fi  =  0,  pure  nested  band)  and  the  6-site  tight  binding  band  [49],  with 
At  =  A0(cosfcT  -  cosfc„),  A0  =  0.1/  and  impurity  strength  V0  =  10t  for  the  former 
band  and  V'o  =  5.3(  for  the  latter. 

The  solutions  of  resonant  frequency  corresponding  to  each  factor  in  Eqn.  3.5 
can  be  tracked  as  a  function  of  separation  R  by  minimizing  2\2  separately.  In 
Fig.  3-1  I  show  the  result  for  the  pure  nested  band  :  the  pentagons  and  open 
circles  identify  different  branches  of  the  splitting  resonance  frequencies;  the  left 
panel  illustrates  the  result  of  R  ||  (110)  inter-impurity  separation  and  the  right 
panel  the  R  ||   (100)  case.  It  is  seen  that  each  factor  Va  corresponds  to  an 
oscillating  function  of  /?,  with  the  factor  determining,  e.g.,  fij,  changing  from  site 
to  site  according  to  whether  the  site  is  even  or  odd.  This  is  due  to  the  strong  R 
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dependence  of  the  components  GQ;  in  the  simplest  case,  R  ||  (110)  and  /i  =  0, 

G3(R,u)  =  G,(R,w)  =  0  but  G0(R,a>)  =  £k  cos{krR/V2)cos(kyR/V^)Gl(uj) 

oscillates  rapidly.  At  R  =  0,  the  problem  reduces  to  the  double-strength  single 

impurity  case;  the  factor  T>i  gives  the  resonant  frequency  fij  (2Vo)  and  the  factor 

Z>2  is  1.  At  large  separation  the  tt{  and  Oj  "envelopes"  are  seen  to  converge  to 

^o  (Vo)  with  a  length  scale  of  a  few  £o  —  10a  for  the  parameters  chosen. 

In  the  R  ||   (100)  case,  the  oscillations  of  the  bound  state  energies  with 
increasing  R  are  not  so  simple,  as  seen  in  Fig.  3-1.  The  one  obvious  simple 
difference  from  the  (110)  case  is  that  the  energy  splittings  vanish  much  faster 
with  distance,  as  expected  from  the  discussion  in  Sec.  3.1.2.  Otherwise  the  short 
distance  behavior  of  the  bound  state  energies  is  complicated.  One  can  check  that 
the  resonant  frequency  closest  to  the  Fermi  level  is  fl%  when  R  =  2  +  4w,  n  integer, 
and  nj"  otherwise. 

In  general,  short  inter-impurity  distance  behavior  is  difficult  to  analyze  ana- 
lytically and  it  is  found  that  in  neither  the  (110)  or  (100)  direction  do  resonances 
appear  at  all  for  R  =  1.  The  hybridization  between  wavefunctions  is  so  strong  in 
these  cases  that  the  picture  of  perturbatively  split  1-impurity  states  breaks  down. 
More  importantly,  the  splittings  are  significant  up  to  quite  large  distances.  Param- 
eters in  Fig.  3-1  are  chosen  such  that  £0  «  10a,  as  seen  from  the  right  panel  of 
Fig.  3-1  where  we  indeed  expect  an  e~fi/&  falloff  according  to  the  previous  section. 
On  the  other  hand,  the  left  panel  of  Fig.  3  1  indicates  strong  interference  out  to 
separations  of  30a  or  more!  This  means  that  if  two  impurities  are  oriented  along 
45  degree  with  each  other,  the  interference  effect  will  be  long  range,  due  to  the 
delocalized  nature  of  the  single  impurity  states  along  gap  nodes. 
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Figure  3  1:  Two-impurity  resonance  energies  Qf2  vs  impurity  orientation. 

3,1,3     The  Interference  of  Two  One-impurity  Wave-functions 

In  the  1-impurity  case,  the  T-matrix  is  given  by  Eqn.  2.26  and  it  is  easy  to 
show  (from  Eqn.  2.3)  that  the  change  of  LDOS  due  to  the  single  impurity  is: 


-«?,",  (r,r,w) 


-WJaae+Jafeffi 


(3.10) 


Quite  generally,  the  Green's  function  can  also  be  expressed  in  terms  of  the  exact 
eigenfunction  i>„(r)  of  the  system  (and  its  conjugate  component  ^*(r))in  with  the 
presence  of  the  impurity  [24,  38] 

i>*n(r)il)n(r) 


j/>*(r)*/)„(r) 


u)  —  Q„ 


(3.11) 


where  the  final  approximation  is  valid  for  a  true  bound  state  with  u  very  close  to  a 
particular  bound  state  energy  Q„,  and  will  be  a  good  approximation  in  the  present 
case  to  the  extent  the  resonances  are  well  defined,  in  the  sense  discussed  previously. 
Comparing  with  the  form  (3.10)  allows  us  to  identify  the  positive  and  negative 
energy  wavefunctions  of  the  single-impurity  resonances  (assuming  Vo  >  0): 


,/>±(r)  =  Z± 


G^(r,w)       u  =  £2(; 


(3.12) 
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where  Z±  are  non-resonant  wave  function  normalization  factors.  Note  that  the 
electron-like  bound  state  eigenfunction  is  directly  related  to  the  off-diagonal  bare 
Green's  function,  while  the  hole-like  wave  function  is  proportional  to  the  diagonal 
bare  Green's  function. 

We  can  follow  the  same  procedure  for  the  two- impurity  Green's  function, 
and  ask  how  the  eigenfunctions  at  a  particular  resonant  energy  are  related  to 
the  single  impurity  wave  functions  we  have  just  found.  Since  the  single-impurity 
resonant  energies  are  different  from  the  2-impurity  energies,  this  aualysis  will  be 
valid  to  the  extent  the  splittings  are  small  compared  to  Qa±.  The  Green's  function 
5G(r,  r)  can  now  be  constructed  from  Eqn.  3.8  and  the  wave  functions  read  off  by 
comparing  with  the  spectral  representation  in  the  same  way  as  in  the  1-impurity 
case.  By  examining  Eqn.  3.8,  it  may  be  shown  that,  depending  on  whether  V,  or 
T>2  is  resonant,  the  wave  functions  thus  extracted  will  be  of  definite  spatial  parity, 
V'„(r)  =  ±4>„{-r).  We  find 
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where  G0{aJ>)  =  G°(r  -  R/2)  ±  G°(r  +  R/2),  and  the  Z'£  are  normalization 
coefficients.  These  are  the  two-impurity  odd  (p)  and  even-parity  (s)  resonant  state 
eigenfunctions  expressed  directly  as  linear  combinations  of  the  corresponding  one- 
impurity  eigenfunctions  i>±  given  in  Eqn.  3.12  With  the  aid  of  the  explicit  forms 
of  those  wavefunctions,  we  can  study  how  the  Friedel  oscillations  induced  by  single 
impurity  interfere.  In  general,  impurity  wavefunctions  include  particle  and  hole 
components  because  of  the  anomalous  scattering  processes  with  amplitude  C?i(R). 
However,  in  some  special  cases  where  (7,(R)  =  0,  for  example,  R  ||  (110),  the 
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Figure  3-2:  The  wavefunction  at  resonance  and  corresponding  LDOS.  Impurities 
are  separated  by  R  =  (6, 6) . 

eigenfunction  become  simpler  and  do  not  mix  particle  and  hole  degrees  of  freedom. 

G?iu(r,cj)       ui  =  Qn 


4U*l  =  z± 


Cr[u(r,u))      ui  =  Qi_ 


and 


(3.14) 


(3.15) 


In  order  to  illustrate  to  what  extent  the  approximated  wavefunction  is 
effective,  I  show  the  real  part.  (1st  column)  and  imaginary  part  (2nd  column) 
of  those  wavefunctions  at  energies  equal  to  the  four  resonant  frequencies  (in  this 
separation,  fl-2±/t  =  ±0.0195,  i11±/t  =  ±0.0075)  explicitly  in  Fig.  3-2.  The  density 
probability,  i.e.,  |*2|  (3rd  column)  and  the  actual  LDOS  (4th  cloumn),  defined  as 


G(r,r,w)    =    ^G°(k,w) 

k 

+  J2  e-<k-k'>-rG°(k,  uJ)fkk,^)G°(k'.  w), 


(3.16) 
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Figure  3  3:  Spectra  on  the  nearest  sites  close  to  the  impurities.  See  text  for  details. 

are  plotted  as  well.  It  is  clearly  seen  that  there  is  a  good  agreement  between  the 

spatial  pattern  of  the  |^|2  and  the  exact  LDOS  calculated  at  each  resonant  energy, 

implying  that  near  each  resonant  energy  the  uoiiresonant  contributions  are  quite 

insignificant. 

3.1.4     Local  Density  of  States,  Realistic  Band  and  Standing  Wave  Condition 

The  distribution  of  LDOS  with  different  inter-impurity  orientations  for  pure 
nested  band  is  plotted  in  Fig.  3-4  and  Fig.  3-5.  There  are  several  novel  features 
which  need  to  be  addressed.  First,  It  is  clear  that  some  states  involve  constructive 
and  some  destructive  interference  between  the  1-impurity  wavefunctions  in  different 
regions  of  space,  but  the  spatial  patterns  are,  not  unexpectedly,  considerably 
more  intricate  than  the  "hydrogen  molecule"  type  states  one  might  first  imagine 
would  form,  with  electrons  living  either  directly  between  impurities  or  completely 
expelled  from  this  region.  This  is  of  course  due  to  the  d-wave  character  of  the 
medium  in  which  the  quasiparticles  propagate.  For  example,  the  LDOS  is  zero 
at  the  point  halfway  between  the  two  impurities  for  the  p-wave  states,  but  it  is 
quite  small  in  the  s-wave  states  as  well.  It  is  furthermore  clear  from  the  figure 
that  both  a  (bonding)  and  p  (antibonding)  functions  can  have  either  constructive 
or  destructive  character,  in  the  molecular  sense.  Note  that  the  states  are  shown 
arranged  vertically  according  to  their  eigeneuergies,  but  recall  that  the  ordering  of 
the  s  and  p  (V^  and  V2)  states  changes  according  to  whether  R  is  even  or  odd,  as 
indicated  in  Fig.  3  1. 
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In  addition,  one  can  look  into  the  spectra  on  some  particular  sites  close  to 
the  impurities.  Fig.  3-3  displays  the  spectra  on  several  sites  whose  positions 
are  illustrated  as  the  numbers  in  the  insets  and  the  impurities  positions  are 
Ri   =   (-3.0), R2  =   (3,0)  (left),  R,   =  (-3,-3),.%  =  (3,3)  (right).  It  is 
surprising  to  observe  that  there  are  some  surprisingly  sharp  peaks,  which  always 
happen  at  high  energies  and  are  far  sharper  than  a  single-impurity  resonance  at  the 
same  energy.  This  is  counterintuitive  based  on  our  knowledge  of  the  one-impurity 
problem:  the  T-matrix  denominator  S-t  possess  an  imaginary  part  which  vanishes 
linearly  when  Fermi  surface  is  approached,  hence  a  resonance  with  lower  energy 
naively  has  narrower  width  (longer  lifetime).  This  novel  peculiarity  must  then  stem 
from  the  quantum  interference  effect.  Indeed,  this  can  be  understood  by  examining 
the  structure  of  the  denominator  of  T-matrix,  with  one  impurity  at  origin: 

V  =  det[V„- 'ts  -  G(Hmp'(R  R, w)]  det f(w), 

where  T  is  the  one-impurity  T-matrix  and  G<hm»>)  is  the  Green's  function  with  one 
impurity  at  the  origin.  Thus  sharp  two-impurity  resonances  occur  for  exactly  the 
same  reason  as  in  the  one-impurity  case,  but  because  the  one-impurity  DOS  at  R  is 
nonmonotonic  in  u,  the  resonance  broadening  is  not  necessarily  proportional  to  the 
resonance  energy. 

One  might  attribute  the  unexpected  "trapped"  (weight  primarily  sitting 
between  impurities)  or  "leaking"  (weight  is  populated  outside  of  the  region  between 
impurities  rather  than  inside)  patterns  to  the  special  symmetries  of  the  half-filling 
tight-binding  band  structure,  for  example,  the  pure  nesting  of  the  Fermi  surface.  It 
is  this  nesting  which  by  allowing  a  considerable  phase  space  for  scatterings  with  a 
fixed,  momentum  vector  Q  (  in  general,  Q  =  (it,  jr)),  leads  to  a  static  or  dynamic 
spatial  modulation  in  terms  of  periodic  arrangement  of  particle  density,  i.e.,  charge 
density  wave  (CDW),  spin  density  wave  (SDW)  etc.  Upon  the  study  on  a  realistic 
band,  we  demonstrate  that  it  is  the  commensuration  of  dominating  scattering 
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Figure  3-4:  LDOS  maps  at  resonant  energies  for  R       ||       (110).  Pure  nested  band, 
Vq  =  lOt;  Impurity  separations  are  shown  on  the  top  of  the  graph. 

wave  vectors  at  the  bound-state  energy  rather  than  the  peculiarity  described  above 
that  determines  the  standing-wave  like  modulation,  a  typical  argument  for  wave 
interference. 

Fig.  3-6  illustrates  the  spectra  on  the  adjacent  sites  of  one  of  the  two  impu- 
rities, with  impurity  separation  R  =  (R,  0)  and  R  =  (R,  R)  for  a  more  realistic 
tight-binding  coefficients  fitted  by  Norman  et  al.  [49]  from  ARPES  data, 

t(k)    =    t0  +  2Ji[cos(fcx)  +  cos(fcj,)]  +AtiCOs(kx)cos{ky) 
+2f3[cos(2fcI)  +  cos(2/£„)] 
+2t4[cos(2fcI)  cos(fcy)  +  cos(kx)  cos(2ky)] 
+U5cos(2kx)cos(2ky)  (3.17) 

with  to...  t6  =  0.879, -1,0.275, -0.087, -0.1876, 0.086  and  \t%\  s0.1488eV.  All  the 
curves  have  been  offset  for  optical  clarity.  The  impurity  resonance  dependence  on 
separation  is  transparent:  the  hybridization  is  strong  when  they  are  too  close  (no 
resonance  pattern  for  R=l);  then  the  splitting  shows  up  and  persists  over  a  certain 
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Figure  3-5:  LDOS  maps  at  resonant  energies  for  R  (100).  Pure  nested  band, 

Vo  =  lOt;  Impurity  separations  are  shown  on  the  top  of  the  graph. 

range  of  inter-impurity  distance  for  each  case.  The  sharp  high  energy  peaks  emerge 
sporadically  and  those  along  (100)  direction  occur  more  frequently  than  that  of 
(110)  direction. 

In  wave  mechanics,  when  two  waves  are  travelling  along  a  string  in  opposite 
directions,  a  standing  wave  can  be  established  as  the  position  and  time  depen- 
dence of  the  resulting  wave  could  be  separated(when  the  two  waves  are  180 
out-of-phase  with  each  other  they  cancel,  and  when  they  are  in-phase  with  each 
other  they  add  together).  In  our  case,  the  formation  of  these  trapped  states  is 
surprisingly  explicable1    in  a  similar  fashion.  For  a  real  bound-state,  the  equation 
G^R,  ui)TG°(R,  ui)T  =  1  should  be  satisfied  strictly  for  both  real  and  imaginary 
parts  at  particular  energies.  The  product  G°(R,u))fCJ°(R,u)f  in  Enq.  3.4  could 


1  It  is  a  bit  surprising  in  that  this  native  picture  doesn't  take  into  consideration 
the  ui  dependence  of  the  T-matrix,  or  the  Nambu  structure  of  Green's  functions. 
Moreover,  the  quantitative  justification  of  the  analysis  is  unexpectedly  obtained 
although  there  are  some  fundamental  difference  between  the  one-impurity  and  two- 
impurity  problems  as  discussed  below. 
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Figure  3  6:  LDOS  spectra  for  realistic  band  and  V0  =  5.3ri  on  nearest  neighbor 
site.  Left  panel:  impurities  at  (-R/2,0)  and  (/J/2,0)  (  R  =  {R,0)),  spectra  taken 
atr  =  (ft/2,1).  Right  panel:  impurities  at  (-ft/2, -ft/2)  and  (ft/2,  ft/2)  (  R  = 
(ft,  ft)),  spectra  taken  at  r  =  (ft/2,  ft/2  +  1). 


be  equivalently  written  as 


£]  eiqRG(k,  ^ffwJGfk  +  q,  w)f  (u 

k,q 


(3.18) 


It  was  argued  by  Hoffman  et  al.  [69,  70] ,  that  the  characteristic  wavevectors  in 
the  spatial  Fourier  transform  of  the  LDOS  of  disordered  BSCCO  samples  are 
determined  by  peaks  in  the  joint  density  of  states  (See  Sec.  3.3.1  for  details),  i.e., 
J2k  Im  Gn (k.  w)Im  G\\ (k  +  q, u),  as  illustrated  in  Fig.  3-16  Although  the  above 
argument  was  proposed  for  single  impurity  and  in  Eqn.  3.18  we  have  a  separation 
dependent  phase  factor  e""R  which  is  deleterious  to  this  argument  by  averaging  the 
joint,  density  of  states  with  random  phase  factors,  it  is  still  found  that  the  long- 
lived  two-impurity  bound-states  are  primarily  derived  from  a  a  few  selected  wave 
vectors;  in  order  to  have  standing  waves,  the  phase  accumulated  from  scattering 
off  from  two  impurities,  together  with  the  contribution  from  space  propagation 
should  sum  up  to  nir,  i.e.,  q  •  R  +  2r;u  =  nir,  with  q  being  the  dominating  vectors 
of  scattering  process  and  </o  the  phase  shift  from  one  impurity.  In  the  unitary 
limit,  2rj0  =  7r,  so  the  comuiensurability  requires  qi  ■  R  +  2r/0  =  q2  ■  R  +  2r/0  = 
Mr,  q3  •  R  +  2t/o  =  mit  to  be  satisfied  simultaneously.  It  is  easy  to  check  for  our 
present  band  structure  that  (3,0),  (7,0),  (11,0),  (14,0)  and  (3,3),  (11,11)  roughly 
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Figure  3-7:  Fermi  surface  of  BSCCO-2212  with  constant  energy  surfaces  at 

u;      =       0.04  shown  as  small  filled  ellipses  at  the  nodal  points,  qi.  q2.  q3  are  wave 

vectors  for  which  the  joint  density  of  states  is  large. 


balance  the  equation  and  we  indeed  see  sharp  resonant  states  for  those  specific 
configurations.  Also,  it  is  observed  that  the  sharp  resonant  states  occur  much 
frequently  in  the  (100)  direction  than  in  the  (110)  direction  as  the  criterion  holds. 
We  also  noticed  that  the  splitting  surprisingly  persists  up  even  to  a  relatively 
large  distance  R  =  13  in  (100)  direction,  while  the  exponential  decay  of  wavefunc- 
tions  in  the  gas  model  would  predict  weak  interference  and  negligible  splitting  at 
this  separation.  This  results  suggests  that  even  with  a  relatively  dilute  concen- 
tration of  impurities,  the  spatial  LDOS  around  impurities  should  not  be  identical 
because  of  the  random  distribution  of  impurities  and  the  interference  between  their 
wavefunctions;  however,  STM  experiments  apparently  observed  point-like,  well 
isolated,  and  nearly  identical  impurities  patterns  embedded  in  the  superconducting 
bath.  In  the  following  sections,  1  will  show  how  this  contradiction  can  be  relaxed. 
Finally,  it  is  also  noteworthy  that  while  the  gas  model  or  pure  nested  band  an- 
ticipates four  distinctive  resonant  peaks,  there  are  usually  only  2  visible  resonant 
peaks  which  can  be  observed  for  the  realistic  band. 
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3.2     DOS  of  two  Impurities  and  Many  Impurities  at  u  =  0  :  Local  Resonance  vs. 
Cumulative  Interference 

3.2.1     Perturbative  Prediction:  SCTMA  and  Its  Validity 

The  problem  of  low-energy  d-wavc  quasiparticle  excitations  has  been  treated 
traditionally  with  a  so  called  self  consistent  T-matrix  approximation  (SCTMA) 
[48,  58,  59,  60]  which  partially  sums  up  the  perturbation  series  and  models  the 
impurities  as  strong  point-like  scatters.  As  illustrated  in  Fig.  3-8,  the  self-energy  is 
approximated  as  E(w)  m  mT(w),  i.e.,  the  product  of  the  concentration  of  impurities 
and  the  single  impurity  T-matrix,  by  dropping  all  the  crossed  diagrams.  It  is  then 
clear  that  SCTMA  incorporates  arbitrary  potential  strength  appropriately  but 
neglects  inter-impurity  correlations.  SCTMA  predicts  that  the  impurity  states 
broaden  with  the  increase  of  impurity  concentration  and  finally  evolve  into  a 
subgap  impurity  band.  The  corresponding  roughly  constant  low  energy  region  of 
the  DOS  is  referred  to  as  the  impurity  plateau,  as  illustrated  in  Fig.  3-9.  While 
the  existence  of  excess  quasiparticle  states  (the  impurity  band)  has  been  confirmed 
by  several  transport  and  thermodynamic  measurements,  the  lifetime  of  those  nodal 
quasiparticles  is  not  consistent  with  what  the  simplest  SCTMA  conjectures.  For 
example,  the  absorptive  part  of  the  microwave  conductivity  a  is  significantly  larger 
than  the  disorder- independent  "universal  limit"  proposed  by  Lee  [71].  Moreover, 
the  expected  T'^-dependence  of  a  at  low  temperatures  (one  factor  of  T  comes 
from  the  density  of  states  and  another  one  comes  from  the  self-energy  which  is 
inversely  proportional  to  T)  [72]  are  not  borne  out  by  experiments  as  well;  instead, 
a  displays  a  linear  T  dependence  in  YBa2Cu307_i  [73,  74,  75]. 

Among  the  various  explanations  proposed  for  these  discrepancies,  one  crit- 
icism casts  doubts  on  the  validity  of  the  SCTMA  itself  by  arguing  that  crossed 
interference  diagrams  neglected  in  this  approximation  leads  to  a  second-order 
n?w2log2aj  divergence  in  2D  coming  from  the  gap  nodes  [64],  to  be  compared  with 
the  SCTMA,  where  one  has  2nd  order  contribution  of  n.;uA  Thus  for  any  fixed  nit 
the  crossed  diagrams  dominate  at  sufficiently  low  energy. 
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Figure  3  8:  Schematic  plot  of  the  self-energy  diagram  with  many  impurities, 
(a)  Zero-energy        (b) 


Impurity  band 

Figure  3-9:  Schematic  figure  of  the  many- impurity  DOS  (a)  in  the  unitary  limit  of 
the  half-filled  band  and  (b)  in  the  unitary  limit  of  a  generic  band. 

3.2.2     Nonperturbative  Predictions: 

When  casting  doubt  on  the  effectiveness  of  SCTMA,  Nersesyan  et  al.  [64] 
accordingly  performed  a  nonperturbative  calculation  for  scattering  within  a  single 
node  and  predicted  a  vanishing  DOS  at  the  Fermi  level,  p(u>)  ~  u°  with  exact 
bosonizatiou.  After  this  work,  several  supposedly  exact  nonperturbative  theories 
made  diverse  conclusions:  p(0)  may  vanish  with  different  power  laws  [61,  76], 
saturate  at  a  finite  value  [65],  or  even  divergefl,  62,  77,  78].  Fig.  3-9  displays 
the  band-dependent  p(uj  =  0)  schematically.  All  these  nonperturbative  works 
claimed  to  yield  exact  results  for  the  problem  of  low-energy  DOS  of  a  d-wave 
superconductor.  These  apparently  conflicting  results  were  reconciled  by  Hirschfeld 
and  Atkinson  [63]  and  Yashenkin  et  al.  [66] ,  who  pointed  out  the  d-wave  disorder 
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problem  was  very  sensitive  to  symmetries  of  both  the  disorder  model  and  the 
normal  state  band. 

Among  all  the  theories,  the  one  proposed  by  Pepin  and  Lee  [1]  is  of  great 
interest  to  me,  which  claims  that  the  result  for  a  N-impurity  T-matrix  is  essentially 
n;  times  the  single-impurity  result.  However,  a  closer  inspection  reveals  that  this 
peculiarity  is  tied  to  the  disorder  interference  rather  than  the  single  impurity 
resonance  at  the  Fermi  level.  Hirschfeld  and  Atkinson  [63],  Yashenkin  et  al.  [66], 
and  Atkinson  e.t  al.  [79]  further  pointed  out  that  this  divergence  is  the  consequence 
of  a  global  particle-symmetry  particularly  possessed  by  the  half-filled  tight  binding 
model. 
3.2.3     Divergence  in  Pure  Nested  Band  with  Unitary  Scatters 

Hirschfeld  and  Atkinson  [63]  have  expended  considerable  effort  to  reconcile 
the  disagreement  among  the  nonperturbative  approaches.  However,  including  their 
works,  none  of  the  studies  above  have  addressed  this  problem  from  a  local  point 
of  view.  The  two-impurity  problem  offers  a  possibility  to  study  this  problem,  to 
observe  how  the  divergent  DOS  in  the  perfectly  nested  case  might  arise. 

As  we  know  from  Sec.  3-1,  the  Green's  function  for  two  d-like  impurities  can 
be  written  as 

G(r,r',u))    =    C°(r-r',w) 

f  ^f;0(r  -  R„w)f,J(w)G°(Rj  -  r,^). 

The  2-impurity  T-matrix  is  a  4  x  4  matrix  (in  the  Hilbert  space  expanded  by  site 
and  spin  indices)  satisfies 

-l 

(3.19) 


V0-'t3-G°((),uj)  -G°(R,u) 

-&>(R,u>)        V0-lf3  -  G°(0,uj) 


It  is  technically  expensive  to  evaluate  of  the  local  Green's  functions  G°(R,  u) 
and  we  can  only  approximate  them  under  some  specific  circumstances.  For  the 
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Figure  3-10:  Schematic  plot  of  bipartite  lattice. 

tight-binding  hall-filled  band  and  the  limit  u)  — >  0,  we  obtained  an  asymptotic  form 
(see  Appendix  A): 

G°(0,w)  =  ~m^7b,  (3.20) 

2        uj 
where  a  =  N/(2itvfv&),  N  =  4  is  the  number  of  nodes,  vp  is  the  Fermi  velocity 
and  «a  is  the  anomalous  quasiparticle  velocity  |  VkAk|,  and  the  cutoff  A  is  of  order 
of  Aq.  The  expansion  in  u>  for  r  =  (m,  n)  depends  on  whether  n  and  m  are  odd  or 
even.  For  the  (even,  even)  case,  we  have 

G°(r,w)  -f  (-1)^  [Gg(0,w)+wO»(r)]  f„,  (3.21) 

where  C'0(r)  is  a  real  function  of  r.  We  find  similar  leading-order  expressions  for 
(m,n)  =  (odd,  odd), 

G°(m,n,u)^ujCo{r)fo,  (3.22) 

while  for  (m,n)  =  (odd,  even)  or  (even,  odd), 

G°(m,  n,  w)  -»  C\ (r)f,  +  C3(r)f3,  (3.23) 

where  Ci(r),  and  C-i(r)  are  real  coefficients  as  functions  of  r  as  well.  This  dis- 
tinction between  even  and  odd  sites  accounts  for  the  oscillatory  nature  of  the 
wavefunctions  for  the  special  case  that  the  Fermi  wavevector  is  commensurate  with 
the  lattice. 

A  concept  of  "  bipartite"  lattice  should  be  introduced  here,  which  divides 
the  square  lattice  into  two  interleaved  subsets,  as  illustrated  in  Fig.  3-10,  where 
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Figure  3-11:  Change  in  p(u)  due  to  impurities.  Impurity  separation:  R  =  (2,2). 

red  and  blue  sports  belong  to  different  sub-lattices  respectively.  A  stunning  effect 
happens  when  the  two  impurities  are  in  the  same  sub-lattice,  i.e.,  R  =  (even,  even) 
or  R  =  (odd,  odd).  In  that  case,  we  found  det  T  diverges  as 


detT- 


[2wG0(R)Gg(0, u))]-'2    R  =  (even, even) 
Gg(0,w)-4  R=  (odd,  odd) 


and  the  correction  to  total  density  of  states 
5p(u)  - 


l/[u  ki2(A/ui)\    R  =  (even,  even) 
2/[u\n2(A/u>)}    R  =  (odd,odd) 


(3.24) 


It  is  worthy  to  pinpoint  the  origin  of  the  divergence  since  the  similar  result  was 
also  obtained  by  Pepin  and  Lee  (PL)[1],  which  claims  that  by  averaging  over  all  the 
possible  impurity  separation  using  the  approximate  form  G°(R,  u)  ~  1/R  out  to  a 
cutoff  ~  t/R,  p(w  — •  0)  diverges  as 


pH 


(3.25) 


M[ln2(A/W)  +  (7r/2)2]' 

where  A  is  a  cutoff.  In  our  model,  for  the  (odd,  odd)  case,  the  local  Green's 
functions  dominate  over  the  nonlocal  ones  and  the  physics  of  the  low-energy 
resonance  is  essentially  that  of  two  noninteracting  impurities,  which  is  verified  by 
the  factor  of  2  in  front  of  the  single  impurity  result(  see  Eqn.  3.24).  For  the  (even 
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Figure  3-12:  Change  in  p(w)  due  to  impurities.  Impurity  separation:  R  =  (2, 1). 


even)  case,  the  interference  effect  seems  to  wipe  out  the  spectral  weight  and  move 
it  to  high  energies  but  the  divergence  survives,  as  its  magnitude  reduces  to  half  of 
that  of  two  isolated  resonances.  In  Fig.  3  11,  I  show  the  change  of  DOS  due  to  two 
impurities  with  separation  R  =  (2, 2)  for  different  impurity  strengths.  It  is  observed 
for  intermediate  impurity  strengths,  four  resonant,  peaks  are  clearly  defined  as  the 
consequence  of  splitting  of  single  impurity  resonances.  When  the  unitary  limit  is 
approached,  a  single  divergent  peak  is  emergent.  However,  it  should  be  stressed 
that,  in  both  cases  of  separation  above,  our  results  are  intrinsically  different  from 
PL  result  which  is  inherently  nonlocal. 

When  impurities  belong  to  different  sublattices  and  are  separated  by  a 
distance  with  site  indices  parity  as  R  =(even,  odd).  The  two-impurity  T-matrix 
defined  in  Eqn.3. 19  is: 


-Gg(0,u)f0  C,(R)fi  +  C3(R)f3 

Ci(R)fi  +  Oi(R)*S         -Gg(0,w)f6 

with  D'  =  G((W)2  -  <7i(R)2  -C3(R)2.  It  follows  easily  that  detf  =  D'~2  and  that 

d 


6p(ui  — >  0)  oc 


duj 


("■*) 


0 


(3.26) 


A  similar  result  holds  for  R  =(odd,  even).  Physically,  the  fact  that  dp  vanishes  at 
the  Fermi  level  indicates  that  bound  state  energies  must  always  arise  at  nonzero 
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Figure  3-13:  (a)  DOS  for  K0=100t.  (b)  Scaling  of  the  DOS  with  V0.  (c)  Scaling  of 
the  DOS  with  L.  (d)  Scaling  of  the  DOS  with  V0  and  L=60. 

energies.  Numerical  calculations  of  the  DOS  shown  in  Fig.  3-12  demonstrate  that 

there  is  no  remnant  of  the  single  impurity  u  ->  0  divergence  for  this  orientation, 

and  that  the  resonance  energies  scale  very  little  with  V0.  In  this  case,  it  is  the 

dominance  of  the  nonlocal  terms  which  shifts  the  resonance  to  finite  energy. 

3.2.4     Discussion:  Bridge  Between  the  Two-impurity  and  the  Many-impurity 
Problems 

We  have  been  aiming  to  show  how  the  two-impurity  problem  is  suggestive  of 
the  fully  disordered  system.  To  explore  their  correspondence,  two  basic  issues  need 
to  be  addressed  here:  1)  what  gives  the  divergence  at  u)  =  0  in  the  pure  nested 
tight-binding  model?  2)  how  does  the  impurity  band  evolve  away  from  w=0? 

6-like  Divergence  at  ui  =  0.  In  Fig.  3  -13,  I  show  the  scaling  of  the  integrated 
DOS  as  a  function  of  the  impurity  strength  and  the  system  size.  The  impurity 
concentration  is  fixed  to  be  n,  =  0.1%.  The  nonperturbative  PL  result  from 
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Eqn.  3.25  is  also  plotted  as  dash  lines  for  direct  comparison.  By  setting  A  =  1, 

we  achieved  a  good  fit  the  numerics  at  V0  =  lOOt.  However,  when  the  impurity 
strength  increases,  the  weight  keeps  concentrating  at  the  center  peak  and  saturates 
slowly.  This  finite-size  effect  requires  us  to  work  in  the  limit  where  the  level  spacing 
is  smaller  than  the  center  peak  width.  The  limit  V0  ->  co  is  illustrated  in  Fig.  3  13 
(d).  Generally,  the  peak  becomes  sharper  when  V0  is  bigger.  It  is  indicated  that  the 
divergence  is  actually  a  delta  function,  i.e.,  limv„_00  p(u)  ~  S(u)  since  the  peaks 
scales  as  p(u)  «  VoF(ojV0). 

What  gives  the  divergence  at  zero  energy  and  why  it  is  not  observed  in  any 
other  many-impurity  case  with  a  generic  band?  It  is  pointed  out  by  Yashenkin  et 
al.  [66]  that  the  half-filled  tight-binding  baud  possesses  a  particular  global  nesting 
symmetry  f2G°(k  +  Q)f2  =  G°(k),  which  in  the  unitary  limit,  i.e.,  U  =  oo  and 
H  =  0  (in  this  this  special  band)  gives  additional  poles  with  momenta  Q  =  (n,  it) 
to  the  particle-hole  and  particle-particle  propagators.  This  diffusive  mode  is  gapless 
and  hence  controls  the  singularity  at  zero  energy.  Any  distortion  on  the  band  or 
deviation  from  unitarity  will  destroy  this  nesting  symmetry  and  yield  zero  DOS 
at  the  Fermi  level.  We  have  performed  careful  examinations  on  the  scaling  of  the 
peak  area  with  respect  to  the  system  size  and  extracted  that  the  localization  length 
£o  ~  40a.  When  system  size  is  bigger  than  £0,  the  localization  effect  makes  the 
effective  system  size  finite  and  the  peak  weight  saturates.  This  is  distinctively 
different  from  the  Pepin  and  Lee's  result  in  origin.  The  divergence  of  the  latter 
arises  from  the  cumulative  effects  of  interference  between  a  large  number  of  distant 
impurities. 

In  Fig.  3-14,  I  show  the  real  space  structure  of  those  states  in  the  divergent 
peak.  The  LDOS  from  the  eigenstates  with  energy  \E„\   <   10~5  is  plotted, 
where  impurities  on  different  sublattices  are  indicated  with  open  and  filled  circles 
respectively.  A  remarkable  discovery  is  that  only  part  of  the  impurities  have 
pronounced  structures  while  others  are  weakly  visible.  It  is  further  noticed  that 
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Figure  3  14:  LDOS  for  2%  concentration  of  impurities  and  |jBn|   <    lfr5i  (20  eigen- 
values). 

those  which  "light  up"  or  "dim  out"  are  located  in  the  same  sublattice  A  or 
B  as  denoted  in  the  previous  section.  While  this  finding  can  be  understood  as 
reminiscent  of  the  two-impurity  problem,  it  is  still  puzzling  in  that  every  impurity 
is  expected  to  have  some  close  neighbours  which  belong  to  a  different  sublattice 
and  hence  contribute  to  the  destruction  of  the  zero-energy  peak,  but  this  is  not 
observed  apparently. 

Away  from  u)  =  0:  A  Real  Space  Perspective  for  STM  Experiments.  An 
amazing  pattern  observed  in  STM  is  that  almost  all  impurities,  regardless  of 
their  local  environment,  appear  to  be  "resonant"  all  through  the  impurity  baud 
[10,  43,  44],  Furthermore,  the  number  of  the  Zn  atoms  matches  the  number  of 
resonances,  which  likely  implies  that  there  are  no  impurity  atoms  "missing." 
These  collective  evidences  above  were  exactly  the  impetus  to  the  proposal  of  single 
impurity  quasiparticle  scattering  scenario.  However,  the  disagreement  between  the 
width  of  resonance  peaks  observed  in  experiments  (which  is  of  the  order  of  ~  10 
meV  and  is  surprisingly  as  the  same  order  of  the  impurity  band  width  in  the  dilute 


unitary  limit  -y  ~  v/ni£>A0!)  and  that  of  theoretical  prediction  on  single-impurity 
(7  <1  meV)  questioned  the  preciseness  of  the  single-impurity  scenario.  Fig.  3-15 
displays  the  spatial  distribution  of  LDOS,  with  a  0.2%  concentration  of  impurity, 
a  value  close  to  the  experimental  observations.  In  the  left,  panel,  LDOS  is  plotted 
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at  a  particular  eigenstate  whose  energy  lies  in  the  impurity  band  but  far  enough 
away  from  the  Fermi  level;  in  the  right  panel,  LDOS  is  plotted  by  averaging  over  a 
final  energy  window  within  which  five  eigenstates  are  included.  It  is  interesting  to 
note  that  while  in  the  left  panel,  only  a  portion  of  impurities  "light  up"  (resonant) 
while  others  are  "dark"  (not  resonant)  simultaneously;  the  right  panel  recovers  the 
experimentally  observed  single  impurity  pattern  qualitatively.  This  possibly  implies 
that  the  DOS  plateau  for  10  >  0  is  formed  by  summing  over  many  impurities  with 
"inhomogeneous  broadening"  (different  impurities  become  resonant  and  turn  off 
at  different  energies,  the  basic  fact  is  that  the  resonant  lifetimes  are  different  for 
different  energies,  since  we  have  a  nontrivial  u-dependent  DOS  for  the  pure  system 
with  which  the  resonant  states  overlap).  We  are  then  inspired  to  formulate  such 
a  hypothesis:  the  exact  image  obtained  by  STM  tips  could  be  the  summation  of 
many  N-impurity  eigenstates,  some  of  which  are  resonant  at  a  given  energy  and 
others  are  not,  over  a  finite  energy  window,  and  the  "single"  impurity  pattern  is 
visually  recovered  by  piling  up  the  snapshots  at  different  energies.  Notice  that  even 
though  the  energy  resolution  in  STM  is  as  small  as  lmeV,  there  are  still  numerous 
eigenstates  within  such  a  window  size  for  a  typical  sample  size  (L  ~  500/1)  and 
impurity  concentration  (n,  ~  0.2  —  0.5%)  with  which  the  tips  may  bin  over  to 
produce  the  nearly  isolated,  four-fold  symmetric,  localized  impurity  patterns.2 
3.2.5     Conclusion 

I  would  like  to  finish  the  discussion  on  two-impurity  problem  by  summing  up 
the  physics  of  most  importance  we  learnt  from  this  specific  model: 

•  The  commensurability  of  the  nodal  wave  vectors  and  the  tight-binding  lattice 
classifies  the  two-impurity  configuration  into  two  sublattices:  the  impurity 
pairs  on  the  same  sublattice  contribute  to  a  divergence  in  p(uj).  At  to  =  0, 


2  The  extended  tails  are  blurred  by  the  incommensurability  between  the  lattice 
and  the  wavevectors  of  eigenstates,  and  this  also  helps  to  set  up  the  single  impurity 
pattern.  Apart  from  the  Fermi  level,  this  phenomenon  is  even  more  pronounced. 
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Figure  3-15:  LDOS  for  0.5%  concentration  of  impurities.  Left:  En  =  0.03854;  right: 
averaged  over  five  eigenvalue  in  energy  interval  \E„  —  0.03f|  <  0.02t. 

with  the  global  constructive  interference  by  the  particular  Ta  symmetry  in 
hall-filled  tight  binding  band,  the  DOS  of  the  many-impurity  system  diverges 
and  contains  a  ^-function  form.  However,  I  should  emphasize  that  this  sharp 
peak  at  a;  =  0  is  the  consequence  of  the  special  feature  of  hamiltonians  with 
Ti  symmetry  only,  and  is  not  a  generic  attribute  of  d-wave  superconductors. 
•  Away  from  the  Fermi  level,  the  quantum  interference  between  impurities  still 
exists.  Incommensurate  scattering  wavevectors  allowed  by  the  increase  in  the 
available  kinematics  phase  space  smears  the  nodal  network  and  distorts  the 
resonant  states,  which  are  inhomogeneously  distributed.  However,  the  finite 
energy  window  of  STM  technique,  although  it  is  quite  small  in  laboratory 
terms,  still  introduces  considerable  homogeneous  broadening  to  the  resonant 
states  by  essentially  averaging  over  a  large  number  of  eigenstates  of  the 
macroscopic  system.  The  isolated  impurity  pattern  with  classic  fourfold 
symmetry  is  visually  recovered. 

3.3     Power  Spectrum  for  Many  Impurities 
When  pebbles  are  thrown  into  water,  they  cause  ripples.  Impurities  induce 
"ripples"  in  the  electron  sea  in  the  similar  way.  If  the  images  of  the  impurity 
induced  "ripples"  are  Fourier  transformed  (FT)  into  momentum  space,  they 
generally  select  out  some  defining  wavevectors  from  which  we  can  map  out  the 
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Figure  3-16:  Left  panekthe  constant  energy  contours  with  the  q  vectors;  right 
panel:  the  evolutions  of  q  vectors  as  a  function  of  energy. 

dispersion  of  eigenstates  in  the  pure  system.  In  normal  metals,  they  are  mainly 
Fermi  wavevectors  [36,  37];  in  a  d-wave  superconductor,  those  Fourier  transformed 
DOS  (FTDOS)  patterns  contain  fruther  useful  information  on  the  SC  state. 
3.3.1      Weak  Scattering  Limit:  Octet  Model  and  Kinematics  of  Quasiparticles 

The  Cornell  STM  group  [69,  70]  first  obtained  the  Fourier  transformed 
density  of  states  for  optimal-doped  disordered  BSCCO,  which  is  defined  as  /)(q)  = 
X]re'qrp(r).  They  plotted  its  absolute  value  as  a  function  of  momentum  q  and  bias 
energy  ui  and  observed  some  distinct  peaks.  An  explanation  based  on  quasiparticle 
scattering  identified  those  peaks  as  the  consequence  of  a  disturbance  due  to  atomic 
scale  disorder.  It  was  speculated  [80]  that  the  scattering  processes  which  have  a 
momentum  shift  q  =  k/iTOj  -  kinitiai  connecting  the  regions  with  substantial  density 
of  states  will  dominate  others,  because  there  is  large  available  phase  space  both  in 
the  initial  and  final  states.  In  a  band  theory,  the  quasiparticle  DOS  p(ui)  at  w  is 
proportional  to 

/      |  Vk1  E(k)\dk,  (3.27) 

JE(k) 

and  the  area  with  the  smallest  energy  gradient  contributes  to  p(uj)  mostly.  If  there 
exist  well  defined  quasiparticles  in  the  superconducting  state,  then  certain  segments 
of  the  Fermi  surface  will  be  gapped  out  and  the  constant  energy  contour  (CEC,  the 
surface  of  which  all  the  points  have  the  same  energy  {k  :  £(k)  =  ui0\)  deforms 
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in  the  d-wave  case  into  "bananas",  as  schematically  plotted  in  Fig.  3  16.  The 

quasiparticle  DOS  around  the  tips  of  bananas  is  higher  than  the  other  domains 

and  we  should  be  able  to  observe  a  pronounced  |p(q)|  for  the  q  vectors  connecting 

these  tips.  In  addition,  these  q  vectors  should  disperse  with  STM  sample  bias  u  in 

a  characteristic  way  as  the  quasiparticle  CEC's  change.  This  proposal  is  called  the 

"octet"  model  for  quasiparticle  scattering  since  for  each  binding  energy  there  are 

eight  identical  q  vectors  connecting  the  "banana"  tips  in  the  first  Brillouin  zone. 

We  can  work  in  momentum  space  within  the  T-matrix  formalism  and  obtain 

the  explicit  form  of  p(q)  (after  taking  off  the  homogeneous  component)  as  well: 

p(q)  ~  lm{J2  G(0) (k  +  q,  u )T(k  +  q,  k, u)Gm (k,  u)] .  (3.28) 

k 

Of  course  the  p(q)  derived  from  the  complete  T-matrix  includes  not  only  the 
octet  peaks  but  also  additional  structures  that  may  hinder  the  visibility  of  the 
desired  ones.  However,  in  the  weak  scattering  limit,  i.e.,  T(k',k)  ~  K(k',k),  may 
scattering  of  high  order  be  neglected  and  the  disorder  potential  be  factored  out  as 
well: 

p(q)  ~  Im{V'(q)  ^C(0)(k  +  q,w)G<0)(k,w)},  (3.29) 

k 

where  V(q)  is  the  Fourier  component  of  the  real  space  potential  K(r).  The  term 
J]kG(0)(k  +  q,a;)G(0)(k,Li;)  is  called  the  response  function  (also  named  A  function) 
which  is  only  determined  by  the  kinematics  of  the  pure  system.  The  physical 
interpretation  of  Eqn.  3.29  is  transparent:  the  imaginary  parts  of  the  Green's 
functions,  namely,  density  of  states,  piniUai  and  pfinai  determine  the  probability 
of  scattering  events  jointly.  Furthermore,  the  magnitudes  of  those  q's  display 
two  distinct  dispersions:  while  some  of  them  increase  with  increasing  bias,  others 
decrease  monotonically  (see  Fig.  3-16).  The  comparison  between  the  experimental 
and  theoretical  w-dependence  of  those  q  peaks  is  presented  in  Fig.  3-16,  showing 
good  agreement  with  each  other!  The  theoretical  prediction  is  calculated  from  a 
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standard  BCS  theory  with  quasiparticle  dispersion 


Ev  =  \Jtl  +  Al  (3.30) 


where  Aj<  takes  dx2-vn  form.  It  should  be  noted  that  Eqn.3.3.1  has  eontributions 
also  from  the  real  parts  of  the  Green's  functions  as  well. 

The  result  of  FT-STS  analysis  has  several  implications:  it  suggests  that  opti- 
mally doped  BSCCO  is  qualitatively  well  described  by  conventional  BCS  quasipar- 
ticles.  Moreover,  because  the  unperturbed  Green's  function  contains  the  dispersion 
of  free  quasiparticles  in  a  homogeneous  system,  we  then  can  trace  the  energy 
dependence  of  those  vectors  and  map  out  the  kinematics  of  the  clean  system!  In 
fact,  the  q  vectors  extracted  from  FT-STS  have  been  used  to  construct  the  Fermi 
surface  of  BSCCO  and  a  comparison  with  that  from  angle-resolved  photoemission 
experiments  displayed  a  good  agreement  [70,  81]  which  confirms  our  knowledge 
of  the  electronic  structure  of  BSSCO  and,  probably  more  importantly,  reinforces 
our  confidence  in  both  the  real  space  (STM)  and  momentum  space  (ARPES)  tech- 
niques because  the  matrix  elements  for  these  physical  processes  are  quite  different. 
It  is  interesting  to  note  that  the  "signal"  used  to  determine  the  Ferrni  surface  of 
he  clean  system  via  this  technique  vanishes  in  the  absence  of  disorder.  Finally,  it 
can  also  be  used  to  measure  the  momentum-dependence  of  the  gap  function  A^.  A 
result  performed  by  Hoffman  et  al.  gave:  A^  =  Aa[Acos(20k)  +  B cos(6(?k)],  with 
Ao  =  39.3  meV,  A  =  0.818,  B  =  0.182,  which  reflects  approximately  the  proposed 
dx2_y2  symmetry  of  the  pairing  function  in  BSCCO. 
3.3.2     Power  Spectrum  for  Many  Impurities  in  BiaSraCaCuaOs 

The  "octet"  model  demonstrates  that  effectiveness  of  the  quasiparticle 
scattering  scenario  within  the  conventional  BCS  theory  for  the  superconducting 
state  of  optimal-  or  over-  doped  samples.  The  study  on  the  Fourier  transform  of 
the  Friedel  oscillations  due  to  impurities  further  provides  important  information  on 
the  kinematics  of  pure  systems.  However,  there  are  also  some  distinct  features  in 
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experimental  data  which  is  not  explicable  within  this  simple  model,  for  examples, 
the  resolved  qi,t}7  peaks  are  much  broader  than  expected;  furthermore,  they  have 
roughly  equal  weight  while  the  octet  model  predicts  the  peak  intensity  of  qr  should 
be  much  stronger  that  of  qj . 

Is  the  quantum  interference  between  the  Priedel  oscillations  responsible  for 
these  discrepancies?  In  the  study  of  two-impurity  problem  [54],  I  have  shown  that 
the  interference  pattern  can  persist  up  to  a  relatively  long  separation;  practically, 
concentration  of  impurities  in  samples  generally  yields  an  average  inter-impurity 
distance  over  which  the  impurities  are  still  quite  "visible"  to  each  other.  If  we 
further  recall  1)  PL's  result  of  the  low  energy  DOS  obtained  by  averaging  over  the 
randomly  distributed  impurity  configurations:  the  interference  gives  an  expression 
for  N-impmity  T-matrix  which  is  essentially  n*  times  the  single-impurity  result; 
2)  the  nearly  isolated  impurity  pattern  could  be  the  artifact  of  homogeneous 
broadening  by  the  energy  window  size  of  STM.  We  should  then  doubt,  the  precision 
of  the  "octet"  model  and  ask  the  following  questions:  1)  is  it  appropriate  for  us 
to  a  priori  assume  that  one-impurity  pattern  survives  from  many  impurities  and 
manipulate  the  experiments  data  on  fully-disordered  samples  within  the  frame  of 
this  theory?  2)  to  what  extent  will  its  predictions  such  as  peak  positions  and  peak 
widths  be  modified  qualitatively  and  quantatively  by  the  interference  effect? 

Many  Random  Weak  Scatters.  Capriotti  et  al.  [82]  and  Zhu  et.  al.  [83]  an- 
swered the  first  question  by  studying  the  power  spectrum  of  N  randomly  dis- 
tributed point-like  impurities.  In  these  works,  the  N-impurity  T-matrix  is  expressed 
in  terms  of  the  1-impurity  T-matrix  tt  =  [1  -  K(7°(r  =  O)]'1^  by 

N, 

Ttj    =    itSiJ  +  ^2it[i-Sm,i}Cl0(Ri-'Rm)fni,  (3.31) 

m=l 

where  the  impurity  potential  at  Rj  is  Vt  =  V0f3,  and  f,  are  the  Pauli  matrices.  In 
the  dilute  concentration  limit,  only  the  leading  order  in  f,  of  Eqn  3.31  needs  to  be 
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Figure  3-17:  FTDOS  at  u   =14  meV  for  weak  potential  scatters  (V0  =   0.67<i):  (a) 
for  one  weak  impurity,  with  a  few  important  scattering  wavevectors  indicated;  (b) 
for  0.15%  weak  scatterers.  Cuts  through  the  data  of  (a) (thick  line)  and  (b)(thin 
line)  along  the  (110)  direction  and  scaled  by  1/-//V7  are  plotted  vs.  qx  in  (c),  while 
(d)  shows  the  weak  scattering  response  function  Im  A3(q,u;).  Peaks  at  q     =     0  are 
removed  for  clarity.  In  all  the  figures,  the  x  and  y  axes  are  aligned  with  the  Cu-0 
bonds. 


considered: 


5f)(q,  w) 


~  Y,  UqJImlc^A^q.u;)] 


(3.32) 


where  t„(q)  =  /„  ^  e-"R-  and  AQ(q,u,)  =  £„  [G0(k,w)fQG°(k  +  q,w)]n  . 
In  the  weak  potential  limit,  Eqn.  3.32  reduces  to 


6p(q,ui)  ~  -K(q)Im  A3(q,w)/7r, 


(3.33) 


which  is  also  valid  for  finite  range  V(r). 

It  is  realized  that  both  of  the  limits  afford  a  complete  separation  between 
degrees  of  freedom  associated  with  the  disorder  potential  and  those  of  the  pure 
system.  There  are  several  consequences  of  this  separation:  1)  the  "octet"  peaks 
are  not  broadened  or  shifted  by  disorder  and  it  is  A„,  the  response  function  of  the 
clean  system  that  determines  the  peak  positions  and  the  widths  in  the  FTDOS; 
2)  the  randomness  of  disorder  introduces  noise  which  is  actually  as  large  as  the 
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signal  since  the  disorder  average  gives  (|d~p(q,w)|2)  ~  Ni  and  (\5p(q,u))\')  - 

{\Sp(q, u)\  )    ~  N](Nj  -  1);  consequently,  some  one  impurity  peaks  will  not  be 
robust  in  many-impurity  case  (for  example,  the  weak  q!  peaks  in  Fig.  3-17  (a)  are 
lost  in  Fig.  3-17  (b));  3)  technically,  we  may  also  compare  the  response  function 
and  the  many-impurity  FTDOS  at  different  energies  to  map  out  the  shape  of  the 
potentials  K(q)  if  they  are  spatially  extended.  Actually,  since  there  is  still  lack 
of  the  knowledge  on  the  spatial  form  of  the  disorder  in  BSCCO,  this  proposes  in 
principle  a  method  to  resolve  the  impurity  structure  from  LDOS  measurement 
directly. 

Finally,  it  should  be  pointed  out  that  the  imaginary  part  of  the  single  particle 
self-energy  due  to  the  point-like  scatters  does  not  broaden  the  "octet"  peaks 
but  suppresses  their  magnitude  in  momentum  space  (at  any  energies)  since  it  is 
momentum-isotropic,  i.e.,  <5£/dk  =  0. 

Unitary  Scatters  and  Extended  Weak  Scatters. 

Point-like  impurities  with  either  dilute  concentration  or  weak  strength  prove 
to  be  inadequate  to  address  on  the  relevant  peak  widths  and  peak  intensities  of 
q!  and  q7  vectors  in  experiments.  Probably,  the  former  fails  because  it  neglects 
the  interference  effect  and  the  latter  further  excludes  the  contribution  from  other 
channels  of  the  response  function.  In  Fig.  3-18,  1  show  r0,  n  and  r3  components  of 
AQ  at  several  energies.  It  is  observed  that,  interestingly,  while  A3  resolves  q7,q2,q6, 
A0  resolves  qi,  q5  quite  clearly  (q3  and  q4  are  present  in  both  of  the  two  channels). 
This  difference  between  A0  and  A3  seems  to  necessitate  the  inclusion  of  unitary 
impurities  as  possible  scattering  centers  and  is  supported  plainly  from  the  structure 
of  the  T-matrix  for  point-like  scatters, 


J0(uj)to_  CT3 

-SoM       c2-.9oM 


1  (U>  -   ,2  _  „2,,.,1   +  r2  _n2,,A  (3-34) 


where  g0{u>  — >  0)  ~  u).  In  the  weak  scattering  limit  the  second  term  in  Eqn  3.34 
dominates  and  qi  and  q5  peaks  are  missing.  However,  in  the  unitary  limit,  i.e., 
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Figure  3-18:  The  plot  of  Nambu  component  of  spectral  functions  at  several  fre 
quencies. 
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C  =  g(u)  — »  0)  -»  0,  r0  and  r3  channels  will  have  roughly  equal  weight  and  therefore 
bring  up  all  the  q  peaks  at  lower  energies.  Additionally,  the  presence  of  unitary 
scatters  is  also  consistent  with  the  finding  of  the  plainly  evident  local  resonances  in 
the  LDOS  from  STM  for  w  <  15meV.3 

Another  constraint  on  the  impurity  potential  imposed  by  the  experiments 
is  the  peak  widths  of  those  vectors.  If  we  allow  that  experiments  may  be  seeing 
background  features  as  well  as  (or  instead  of)  the  true  "octet"  peaks,  then  the 
FTDOS  in  Fig.  3-17  is  fairly  compatible  with  experimental  data  at  -14  meV  (as 
shown  in  Fig.  3-19).  However,  comparisons  at  other  energies  are  significantly 
worse.  As  revealed  by  the  experiments,  besides  the  observance  of  unitary  scatters, 
the  energy  gap  extracted  from  STM  measurements  of  BSCCO  is  inhomogeneous 
at  the  nanometer  scales  [84,  85],  ranging  from  underdoping  to  slightly  overdoped. 
With  the  lack  of  understanding  of  the  source  of  those  inhomogeneities,  we  modified 
them  as  smooth  random  potentials,  probably  originating  from  charge  inhomo- 
geneities from  BiO  layer,  with  the  structure  as  V(r)  =  £V ,  V(i)exp(-fj/A)/f;  and 
fj  =  [(r  -  R,)2  +  (fz]i/2,  where  R*  +  Mz  are  the  defect  locations,  V(i)  are  the  defect 
potentials  and  A  is  a  screening  length.  These  extended  impurities  do  not  induce 
any  resonance  feature  at  all  as  they  never  make  real  poles  for  the  denominator  of 
T-matrix,  but  their  momentum-dependent  potentials  V(q)  assign  a  considerable 
probability  to  the  "off-shell"  scattering  (the  process  involved  with  the  k's  which  do 
not  satisfy  the  poles  equations  for  the  Green's  functions  of  initial  and  final  states 
simultaneously)  which  act  to  blur  the  octet  peaks. 


3  This  is  also  true  for  magnetic  impurities  since  their  spin  resolved  matrix  is 
written  as  r0  in  Nambu  notation.  But  we  focus  on  the  potential  scatters  here. 
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In  Fig.  3-19,  I  show  the  comparison  between  theoretically  calculated  FT-DOS 
and  experimental  results.  The  theoretical  results  were  obtained  by  solving  the  Bo- 
goliugov  De-Gennes  (BDG)  equations  ''   including  both  unitary  and  weak,  extended 
scatters.  The  q  vectors  from  the  "octet"  model  are  also  plotted.  A  fair  agreement 
between  the  calculated  and  experimental  FT-DOS  is  achieved  at  low  energies.  The 
qi  and  qs  peaks  are  well-resolved  theoretically,  but  most  of  the  structure  comes 
from  a  set  of  broader  background  features  which  disperse  along  with  the  octet  vec- 
tors simultaneously.  In  general,  the  background  structures  disperse  qualitatively  as 
one  expects  in  single  impurity  model.  Meanwhile,  theoretically  calculated  FTDOS 
display  an  asymmetry  between  positive  and  negative  energies  where  experimental 
ones  do  not.  This  is  maybe  the  result  of  the  large  asymmetry  of  the  model  band,  in 
which  there  is  a  Van  hove  singularity  at  uj  «  50  meV  that  is  not  observed  in  STM 
experiments. 

While  qi .  q-  arc  identified  as  "remnant"  of  octet  peaks  together  with  the 
dispersing  back  ground,  q7  seems  to  be  misidentified  as  octet  peaks  in  experimental 
data,  as  we  argued,  since  the  measured  magnitude  is  as  twice  as  the  predicted 
value  from  nearest-neighbor  rf-wave  model.  Because  the  q7  peaks  comes  from 
intranodal  scattering,  they  are  a  direct  measure  of  the  k-dependenee  of  the 
superconducting  gap  and  scales  with  ~  l/e&,  where  v&  is  the  gap  velocity  at  the 
nodes.  McElroy  et  al.  [70]  have  exploited  this  mechanism  to  map  out  the  gap 
function  At  and  were  forced  to  introduced  a  significant  subleadiug  cos  60  harmonic 
to  fit  their  data.  This  is  in  sharp  contrast  to  the  pure  near-neighbor  d-wave  form  at 
optimal  doping  resolved  from  ARPES  spectroscopy.  Recognizing  that  the  observed 


4  In  BDG  formalism,  the  hamiltonian  for  a  superconductor  on  a  square  lattice  is, 

where  the  angle  brackets  indicate  that  site  indices  i  and  jare  nearest  neighbors,  C/,- 
is  the  impurity  potential,  and  <Vj  =  -V (c^c^)  is  the  mean-field  order  parameter. 
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Figure  3-19:  Comparison  of  theory  and  experiment.  Left  panel:  LDOS;  middle 
panel:  calculated  FTDOS;  right  panel:  FTDOS  from  STM  experiments. 


feature  at  roughly  twice  the  true  q7  is  in  fact  the  background  feature  found  in  our 
calculations  may  enable  one  to  bring  the  two  experiments  in  closer  agreement. 

At  energies  \u.<\  >  15  meV,  Fig.  3  19  shows  a  qualitatively  similar  nanoscale 
variation  of  LDOS  as  seen  in  experiments.  However,  the  (110)  peaks  associated 
with  forward  scattering  are  stronger  than  (100)  peaks.  This  can  be  understood 
since  the  unitary  scatters  are  only  resonant  at  low  energies  and  do  not  have 
noticeable  effects  at  high  energies.  Therefore,  the  difficulty  to  bring  up  qi  peaks 
places  further  constraint  on  the  impurity  potential.  Alternatively,  it  has  also  been 
proposed  that  qi  peaks  could  be  attributed  either  to  the  exotic  orders  such  as 
stripe  formation  [30,  31]  or  to  the  scattering  from  gap  inhomogeneities,  i.e.,  the 
contribution  from  T|  channel. 


CHAPTER  4 

DISCUSSION:  QUASIPARTICLES,  DISORDERS  AND  THEIR  IMPLICATIONS 

FOR  STM  SPECTROSCOPY 

The  notable  improvement  in  the  resolution  of  STM  experiments  has  shifted  the 
focus  of  the  study  of  high-temperature  superconductivity  from  bulk,  macroscopic 
property  to  nanoscale,  specifically,  local  quantum  states  and  electronic  structures, 
a  level  from  which  we  believe  the  fundamental  mechanism  of  high  temperature 
superconductivity  may  eventually  be  uncovered.  The  illuminating  results  of  STM 
experiments  have  established  several  undisputed  facts:  a)  the  spatial  modula- 
tions of  the  LDOS  and  order  parameters  are  inhomogeneous  at  the  nanoscale  in 
Bi2Sr2CaCu-208  ;  b)  disorder  with  unknown  origin  and  structure  is  ubiquitous  in 
this  sample;  c)  optimal-  and  over-doped  cuprates  may  be  well  described  by  the 
conventional  BCS  theory  with  unusual  order  parameters.  I  attempt  to  elucidate 
to  what  extent  we  can  relate  such  atomic-scale  electronic  phenomena  to  disorder 
induced  wave-like  quasiparticle  states  and  the  quantum  interference  effects  be- 
tween them.  I  here  briefly  summarize  my  understanding  on  the  implications  of  the 
quasiparticle  scattering  model  for  the  recent  discoveries  in  STM  experiments. 

4.1      Atomic  Level  Resonances,  Quaspartiele  States  and  Unitary  Scatters 

The  zero-bias  differential  conductance  map  on  the  surface  of  BSCCO  displays 
localized  resonant  features  of  size  of  3  nanometers  and  four-fold  symmetric  distribu- 
tion of  LDOS.  In  addition,  the  spectra  on  the  centers  of  those  bright  spots  exhibit 
a  zero-bias  sharp  peak,  as  shown  in  Fig.  4-1  [43]. 

There  has  been  great  deal  of  effort,  ascribing  these  observations  to  the  quasi- 
particle scattering  from  atomic-scale  defects.  While  the  identities  of  these  "native" 
defects  are  still  not  clear  (  the  possible  candidates  include  defects  such  as  Cu  and 
O  vacancies,  crystal  defects  etc),  the  spectroscopic  signature  are  quite  similar  to 
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Figure  4-1:  Zn  impurity  resonance,  left  channel:  dl/dV  Vs.w;  right  channel:  the 
LDOS  spectrum  above  the  Zn  site. 


those  observed  when  Zn  and  Ni  are  explicitly  substituted  for  Cu.  Theoretical  cal- 
culations, in  which  Zn  and  Ni  were  modelled  as  potential  and  magnetic  scatterers 
respectively,  seemed  to  support  this  conjecture  by  yielding  qualitative  agreement 
with  experimental  results  and  further  substantiated  the  possibility  of  using  local- 
ized defects  as  atomic  probes  for  the  superconducting  state.  The  validity  of  this 
scenario  was  afterwards  reinforced  by  its  successful  application  to  the  analysis  of 
Fourier  transformed  LDOS,  through  which  we  mapped  out.  the  kinematics  of  the 
pure  superconducting  state  (such  as  the  Fermi  surface,  the  band  dispersion  and  the 
momentum-dependence  of  the  order  parameter)  that  was  confirmed  by  the  results 
from  momentum  sensitive  techniques  as  well.  Additionally,  transport  measurements 
[86,  87]  also  demonstrated  the  existence  of  excess  low-energy  quasiparticle  states, 
which  may  result  from  the  collective  interference  between  impurity  wavef unctions. 
However,  there  remain  a  few  difficulties  in  the  quasiparticle  scattering  scenario: 
•  The  expected  spectral  weight  distributions  around  Zn  impurities  are  com- 
pletely reversed  in  experiments  (  see  Chapter  2  for  details).  A  proposal 
emphasizing  the  "filter"  effect  arising  from  the  coupling  between  the  orbitals 
with  zero  in  plane  components  and  the  3dr2^y2  wavefunctions  of  the  nearest 
neighbor  Cu  reconciles  this  disagreement  but  raises  another  debate  on  the 
sign  of  the  Zn  impurity  itself.  So  far,  hand-waving  arguments  based  on  the 
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atomic  configuration  suggests  Zn  is  an  attractive  potential  for  electrons 
while  calculations  determining  the  unitary  limit  by  adopting  realistic  band 
structure  and  ab  initio  DFT  calculations  find  a  repulsive  potential.  This, 
together  with  the  fact  that  Zn  (presumably  a  potential  scatterer)  atoms  are 
more  destructive  to  local  coherence  than  Ni  (presumably  a  magnetic  scatter), 
remains  an  unsettled  issue. 
•  The  ability  of  one-impurity  quasiparticlc  scattering  model  of  any  kind 
to  explain  the  details  of  local  STM  spectra  in  samples  with  percent  level 
disorder  is  severely  limited.  Despite  the  fact  that  conventional  scattering 
theory  might,  be  intrinsically  deficient  since  it  is  completely  blind  to  the 
possible  strong  correlations  between  the  impurity  states  and  the  surrounding 
electrons,  it  is  inadequate  to  address  several  experimental  observations. 
For  example,  although  the  analysis  of  Fourier  transformed  LDOS  (within 
the  single  impurity  scattering  scenario)  resolved  the  characteristic  peaks 
associated  with  the  kinematics  of  pure  system,  the  actual  heights  and  widths 
of  these  peaks  are  inconsistent  with  the  one-impmity  model.  With  the 
inclusion  of  extended  impurities,  proposals  based  on  the  quantum  interference 
between  impurities  yielded  better  agreement  with  the  experimental  data  but 
still  could  not  reproduce  all  the  observations  consistently. 
4,2      Nanoscale  Gap  Inhomogeneity 
In  addition  to  aforementioned  zero  bias  resonances  due  to  planar  defects,  a 
further  fascinating  aspect  of  the  local  electronic  structure  is  the  observation  of 
large,  nanoscale  fluctuations  of  the  LDOS.  In  Fig.  4-2,  I  show  the  plot  of  gapmap 
in  an  overdoped  BSCCO-2212.  In  this  plot,  each  spot  is  a  real  lattice  point  and  the 
color  represents  the  relative  magnitude  of  the  local  gap  maxima  which  is  defined 
as  the  position  of  coherent  peak  in  the  local  spectral  weight  on  the  chosen  site. 
Generally,  this  nanoscale  inhomogeneity  appears  at  biases  close  to  the  bulk  gap 
edge  in  BSCCO  accompanied  by  homogeneous  low  energy  (\ui\  <  A0)  spectra  across 
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Figure  4-2:  The  gapmap  of  overdoped  BSSCO.  Sample  size  500Ax500A. 
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Figure  4-3:  Spectra  on  different  sites  along  the  horizontal  cut. 
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the  sample  (See  Fig.  4-3).  This  observance  was  interpreted  as  the  local  fluctuation 
of  the  order  parameter  and  divided  the  system  into  "patches"  with  fluctuating 
phase  coherence,  which  motivated  the  study  of  "granular  superconductivity."   ' 
However,  we  should  keep  in  mind  that  the  STM  tips  measure  quasiparticle  excita- 
tions directly  rather  than  the  order  parameter;  additionally,  in  a  fully-disordered 
superconductors,  quantum  interference  effects  can  also  give  birth  to  the  local  gap 
inhomogeneity.  While  we  should  question  whether  this  inhomogeneity  is  intrinsic  in 
BSCCO,  we  may  also  ask  how  for  can  we  proceed  with  quaiparticles  and  disorder 
alone,  i.e.,  is  this  phenomenon  completely  explicable  in  "QP+disordcr"  scenario,  or 
is  it  some  intrinsically  exotic  state  that  can  be  distinguished  from  a  highly  disor- 
dered superconductor?  There  was  an  attempt  [29]  to  study  the  second  possibility 
by  using  Ni  to  induce  local  resonance  in  superconducting  regions  and  discriminate 
them  from  regions  of  other  identities.  It  was  found  that  Ni  atoms  only  "light  up" 
in  the  region  with  relatively  small  gap  value  and  higher  coherent  peak  (i.e.,  good 
superconducting  grains)  but  how  the  Ni  resonances  correlated  with  the  inhomo- 
geneity was  not  resolved  explicitly  yet.  As  far  as  the  first  possibility  is  concerned, 
the  homogeneous  shape  of  the  low  energy  spectra  excludes  a  large  concentration  of 
short-range,  strong  scatters  since  they  should  scatter  the  nodal  quasiparticles  and 
result  in  excess  weight  around  the  Fermi  level  presumabl;  furthermore,  the  fact  that 
these  inhomogeneities  take  place  at  energies  close  to  the  bulk  gap  edge  suggests 
that  it  is  the  antinodal  quasiparticles  that  suffer  intense  scattering  from  weak, 
smooth  potentials.  This  could  be  particularly  true  in  the  optimal-  to  overdoped 
BSCCO  where  we  have  a  better  global  superconducting  phase  coherence.  In  Fig. 
4-4  and  Fig.  4-5,  I  show  the  simulation  on  the  gapmap  and  the  spectra  along  a 
horizontal  cut  by  solving  the  Bogoliubov  De-Gennes  equations. 


1  The  superconducting  patches  are  not  granular  in  a  structural  sense.  We  call 
them  "granular"  to  refer  the  fluctuation  of  the  order  parameters  which  characterize 
the  phase  condensate  in  each  "grain." 
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Figure  4-4:  Gap  map  from  numerical  solution  of  BDG  equations  with  only  smooth 
disorder. 


Figure  4-5:  LDOS  along  the  horizontal  cut  in  the  top  figure. 

Finally,  there  is  another  piece  of  evidence  supporting  the  significant  scattering 
of  antinodal  quasiparticles  regardless  of  the  origin  of  the  disorder  [70]:  the  con- 
ductance map  shows  that  for  a  given  uj,  the  observed  Umklapp  scattering  induced 
LDOS  modulation  is  extremely  localized  to  the  region  whose  local  gap  value  is 
approximately  equal  to  u.  This  indicates  that  for  k  ss  (jr,  0),  the  quasiparticle 
dispersion  experiences  strong  nanoscale  variation  and  the  antinodal  quasipaticle 
lifetime  is  remarkably  influenced  by  the  scattering.  Nevertheless,  this  is  in  con- 
tradiction with  the  fact  that  the  antinodal  quasiparticles  are  well-defined  in  the 
superconducting  state,  indicated  by  the  the  sharp  (it,  0)  quasiparticle  spectra  in 
ARPES. 
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4,3     Long-range  Modulation  and  Exotic  Orders 

The  electronic  structure  of  cuprates  has  been  proven  to  be  extremely  sensitive 
to  transitions  into  a  variety  of  ordered  states.  Doping  increases  the  mobility  of 
the  itinerant  holes  and  they  redistribute  among  the  disordered  potential  landscape 
to  establish  significant  electronic  modulations.  Various  spatial  modulations,  such 
as  inhomogeneous  gap  distribution,  incommensurate  periodicity  oriented  along 
the  copper-oxygen  bond,  "checkboard"  pattern  were  consecutively  reported  in 
various  underdoped  samples  [88,  89].  Identification  of  those  spatial  modulations 
generally  falls  into  two  categories:  it  could  be  the  ordering  either  tied  to  the 
lattice  or  tied  to  the  Fermi  surface  (interference  effect).  The  first  possibility  tries 
to  attribute  the  experimental  observations  as  the  consequence  of  competition 
between  the  superconducting  state  and  exotic  orders  such  as  cf-density-wave 
[18],  antiferromagnetic  spin  density-wave  [90]  and  centers  the  debate  onto  the 
exploration  on  the  underlying  mechanism  of  the  "pseudogap"  state,  the  region 
where  the  observations  described  above  are  prominent.  Several  authors  have 
studied  the  impurity  effect  in  those  hypothesized  orders  (or  a  mixed  states  of  the 
exotic  orders  and  the  superconducting  state)  in  the  "pseudogap"  state  [91,  92,  93] 
and  proposed  that  the  spectroscopic  feature,  such  as  the  spatial  distribution 
and  energy  of  the  resonant  states  could  be  used  to  identify  the  different  types 
of  orders.  The  second  scenario,  although  depending  on  the  detailed  information 
about  the  quasiparticle  excitations  for  the  proposed  order,  generally  predicts 
energy-dependent  characteristic  scattering  vectors,  which  is  in  conflict  with  the 
energy-independent  wavevectors  extracted  from  the  Fourier  transformed  STM  DOS 
at  q  =  4.5u0  in  underdoped  samples.  While  it  might  be  premature  to  reject  the 
quasiparticle  interference  scenario  entirely,  this  contradiction  at  least  imposes  a 
constraint  on  the  Green's  function  by  asking  for  a  dispersive  imaginary  part  and 
a  non-dispersive  convolution  between  two  imaginary  components  (similar  to  the 
argument  of  the  joint  density  of  states  )  simultaneously.  As  for  now,  the  content 
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in  this  section  still  remains  the  most  intriguing  puzzle  in  the  understanding  of  the 
phase  diagram  of  cuprates. 


CHAPTER  5 

ANGLE-RESOLVED  PHOTO-EMISSION  SPECTRA  OF  DISORDERED 

CUPRATES 

In  the  previous  chapters,  I  have  discussed  how  the  recent  improvements  in 

STM  technique  facilitated  our  comprehension  of  the  cuprates  by  providing  detailed 

real  space  image  of  their  surfaces  and  extracting  local,  subatomic  electronic 

structures  of  superconducting  materials  directly.  However,  spectroscopies  in 

momentum  space  provide  complementary  information  and  are  sometimes  even 

more  important  because  they  interpret  the  dynamics  and  kinematics  of  the 

systems  in  a  more  "physically  transparent"  way.  The  angle-resolved  photoemission 

(ARPES)  technique  is  such  one  desirable  method  since  the  coherent  part  of  its 

momentum-resolved  spectrum  elucidates  the  dispersion  and  self-energy  of  single 

particle  unambiguously  hi  quasi  two  dimensional  systems.  ARPES  has  been 

empirically  useful  in  determining  the  electronic  states  of  cuprates,  for  example, 

the  energy  distribution  curves  (EDC)  of  ARPES  follow  the  hypothesized  rf-wave 

like  momentum-dependence  of  the  order  parameter  Ak.  Recently,  the  dramatic 

improvement  in  the  angular  and  momentum  resolution  of  ARPES,  as  of  2  meV 

and  0.2°  respectively,  makes  this  technique  a  leading  tool  in  the  study  of  high-T,. 

superconductors.  The  one  particle  spectral  function  resolved  from  the  advanced 

ARPES  experiments  reveals  novel  and  perplexing  line  shapes  [33,  94,  95,  96,  34] 

in  the  low  energy  regime  and  each  has  ushered  in  a  new  round  of  investigation  in 

high-Tc  superconductors.  I  will  concentrate,  in  this  chapter,  on  the  impact  of  elastic 

forward  scattering  and  its  implications  on  the  contemporary  dispute,  regarding  the 

energy-  and  momentum-dependence  of  the  single  particle  self-energy  [97,  98). 


7(> 


77 
5,1     Angle-resolved  Photoeinission  Spectra  of  High  Temperature  Superconductors 
5.1.1     Spectral  Function  and  Self-energy 

When  photons  hit  metal  surface,  they  will  be  absorbed  by  electrons  and  hence 
pump  the  electrons  up  into  higher  energy  levels.  The  excited  electrons  will  then 
propagate  in  the  metals  and  eventually  be  emitted  through  the  surfaces  if  they 
have  sufficient  energy  to  overcome  the  work  function.  The  excitation  of  electrons 
is  governed  by  the  Fermi's  golden  rule.  The  intensity  of  the  photoelectrons  are 
measured  by  the  detectors.  If  we  assume  that  the  relaxation  time  of  the  excited 
electrons  is  much  longer  than  the  time  for  them  to  reach  the  surface  (this  is  the  so 
called  sudden  approximation),  we  could  roughly  formulate  the  intensity  of  outgoing 
electrons  beams  as: 

l(k,u>)  =  I0(k)f(v)AQsi,w),  (5.1) 

where  /0(k)  is  determined  by  the  momentum  dependent  matrix  elements  and  the 
Fermi  function  J(ui)  illustrates  that,  this  process  can  only  measure  the  unoccupied 
states.  The  one  particle  spectral  function  A(k,  ui)  describes  the  probability  of 
creating  or  annihilating  one  particle  in  an  interacting  many-body  system  and  is 
associated  with  the  imaginary  part  of  the  retarded  Green's  function,  as  defined  in 
Eqn.  5.2: 

A(k,tv)    =    ^|<+1|cJ|*f|2*(w-^[+1  +  £f) 

m 

+  Ei^'ki^i^-^-'  +  O 

m 

=    — ImG(k,w  +  i04).  (5.2) 

7r 

In  a  non-interacting  system,  A(k,  uj)  is  nothing  but  a  delta  function,  <5(lj  -  tk), 
everywhere,  due  to  the  infinite  lifetime  of  electrons  in  the  absence  of  inter-body 
interactions.  When  the  many-body  interaction  is  turned  on  adiabatically,  the  bare 
electrons  will  be  "dressed"  by  the  medium  and  can  be  described  by  the  Landau 
Fermi  liquid  theory  (FL)  provided  they  still  live  long  enough  to  be  considered 
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as  approximately  independent  objects.  The  FL  theory  painted  the  big  picture  of 

interacting-electron  dynamics  and  introduced  an  important  theoretical  concept, 
i.e.,  the  "quaisparticles."  With  this  concept,  we  can  treat  the  dressed  "electrons" 
as  bare  entities  with  finite  but  considerable  long  lifetimes.  The  one  particle 
Green's  function  is  generally  modified  as  G^fcw)  =  GjT1  -  E(k,w),  where 
E(k,w)  =  S'(k,w)  +  ffi"(k,w),£'(k,w)  =  Re£(k,w),£"(k,a;)  =  Im£(k,w)  and 
£  is  called  the  "self-energy."  Its  real  part  describes  the  renormalization  of  the  bare 
electronic  dispersion  due  to  many-body  effects  and  its  imaginary  part  is  related  to 
the  one  particle  lifetime  by  r(k,ui)  =  -l/2ImE(k,u>).  The  general  expression  for 
the  spectral  function  is  then: 

\n      \       1  E"(k,w) 

A(k,  u)  = —, i-J — l /cos 

7r[o;-ek-E'(k,u;)]2  +  [E"(k,a;)]2  [    J) 

The  self  energy  £(k,w)  arises  from  many-body  and  impurity  effects  and 
encapsulates  aU  the  information  about  interactions.  Generally,  we  can  plot  the 
measured  intensity  (oc  A{k,Lo))  either  as  a  function  of  u  (binding  energy)  for  a 
fixed  momentum  k  (EDC),  or  as  function  of  k  with  a  fixed  value  of  u  (MDC).  The 
detailed  forms  of  self-energies  extracted  for  the  cuprates  from  these  two  methods 
exhibit  quite  unusual  features  which  are  in  conflict  with  the  traditional  FL  theory 
and  will  be  the  focus  point  in  the  discussions  below. 
5.1.2     Self-energy  in  Normal  and  Superconducting  States 

The  spectral  properties  of  single  particles  in  the  normal  states  obtained  from 
ARPES  data  are  extremely  unusual.  Despite  the  existence  of  well-defined  Fermi 
surfaces  in  various  samples  [34],  the  energy  and  momentum  dependence  of  the 
spectral  function  deviates  considerably  from  the  predictions  of  the  FL  theory. 
The  most  striking  discovery  is  that  there  exist  no  well-defined  quasiparticles  in 
the  normal  state  of  optimally  doped  or  slightly  overdoped  HTS's  near  the  (ir.O) 
point  of  the  Brillouin  zones  [99,  100,  101],  while  in  superconducting  states  of 
these  samples  quasi-particles  are  well-defined  over  a  certain  range  in  the  1st  BZ, 
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Figure  5-1:  ARPES  spectra  from  overdoped  Bi2212  (Tf=87K).  a)  normal-  and  b), 
(c)  superconducting  state  data  measured  at  the  k  point  indicated  in  the  iuset  of 
(c). 

as  shown  in  Fig.  5  1  [100].  The  photoemission  intensity  exhibits  (panel  (b))  a 
pronounced  peak,  followed  immediately  by  a  "dip"  and  then  a  "hump"  feature 
at  higher  binding  energies  in  the  superconducting  state1  .  When  the  temperature 
increases,  the  coherent  low  energy  quasipartiele  peak  evolutes  continuously  and 
finally  melts  into  a  quite  broaden  background  (panel  (a)),  indicating  that  a 
considerable  self-energy  develops  across  the  phase  transition  and  the  quasiparticles 
are  no  longer  well-defined. 

We  might  cast  a  doubt  on  the  applicability  of  the  FL  theory  to  the  normal 
state  of  high  Tc  superconductors  since  the  normal  state  temperature  is  of  the  order 
of  100K  (  it  is  quite  high  compared  to  the  general  definition  of  "low"  temperature), 


1  It  is  called  the  "peak-dip-hump"  feature.  The  high  energy  humps  are  specu- 
lated generally  to  be  the  consequence  of  phonou  modes,  magnetic,  collective  mode 
or  bilayer  splitting  effects. 
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but  we  should  also  not  abandon  it  unthinkingly  because  1)  it  has  assisted  our 
understanding  in  various  crystalline  solids  and  been  useful  over  a  wide  range  of 
temperatures  2)  there  are  still  well-defined  quasiparticles  in  certain  segments  of 
Fermi  surfaces  when  the  gap  opens,  which  is  consistent  with  the  conventional 
BCS  theory.  However,  an  estimate  of  the  full  width-half  maximum  (FWHM)  of 
the  broadened  peak  yields  a  value  of  ~  lOOmeV,  which  is  one  order  of  magnitude 
bigger  than  the  thermal  broadening.  This,  together  with  the  continuous  evolution 
of  the  qp  peak  suggests  that  temperature  cannot  be  responsible  for  such  a  large 
peak  width  directly. 

While  there  are  numerous  scenarios  concerning  the  origin  of  the  one  particle 
self-energy  such  as  electron-phonon  coupling,  spin-fluctuations,  magnetic  resonance 
etc,  a  Marginal  Fermi  liquid  hypothesis  proposed  by  Varma  et  al.  [8]  seems  to 
describe  the  transport  and  thermodynamic  properties  and  the  anomalies  observed 
in  ARPES  phenomenologically.2    It  was  further  speculated  [9]  that  in  order  to 
interpret  the  behavior  of  self-energy,  we  need  include  a  term  independent  of 
frequency  and  temperature,  which  is  probably  associated  with  elastic  forward 
scattering  of  quasiparticles  by  impurities  located  away  from  the  Cu02  plane,  and 
a  MFL  component  as  defined  in  Eqn  5.4.  The  scattering  rate  in  the  normal  state 
then  reads: 

r(k,a>)  =  r„(k,w)  +  Tmpl%w),  (5.5) 


2  The  MFL  theory  was  postulated  to  reconcile  the  conflict  between  the  anoma- 
lous (but  universal)  properties  in  the  normal  state  of  cuprates  and  the  well-defined 
Fermi  surface  (in  optimally  doped  samples).  It  assumes  that  electrons  couple  with 
a  bosonic  excitation  whose  spectrum  has  a  scale-invariant  form  B(uj)    ~   min(j,,  1). 
The  electron  self  energy  due  to  the  scattering  from  this  bosonic  spectrum  has  a 
form: 

£(k,  u)  =  \[ui  log i-x],  (5.4) 

wc        2  ' 

where  x  =  max{\u\,  T).  This  singular  behavior  of  self-energy  leads  to  the  absence  of 
the  FL-like  quasiparticles  even  on  the  Fermi  surface  when  T=0. 
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The  possible  applicability  of  elastic  forward  scattering  is  further  reinforced  by  the 

following  fact:  in  Born  approximation,  the  self-energy  due  to  impurity  scattering  is 
written  as 

E(k,w)  =  n,]T|Vk,kf  C0(k',uj).  (5.6) 

k' 

Qualitatively,  if  V'kk.  is  peaked  at  certain  direction,  then  X(k,u)  is  roughly 
proportional  to  the  imaginary  part  of  the  integrated  Green's  function  over  a 
narrow  range  centered  around  k,  i.e.,  G{k,u)  ~  Ek'-kiJt^C'''^)-  wllen  W  is 
extremely  small,  namely,  in  forward  limit,  E(k,w)  is  proportional  to  !'p(k)_1  (the 
inverse  of  group  velocity  of  band  electrons),  which  generally  takes  the  maximum 
value  at  antinodes  and  minimum  value  at  the  nodes  on  the  Fermi  surface.  This 
momentum  anisotropy  of  E(k,  10)  coincides  with  the  experimentally  extracted 
elastic  component  Fel(k,  uj)  by  assuming  the  forward  scattering  scenario,  which 
also  displays  the  similar  behavior  with  a  maximum  value  of  120  meV  at  antinodal 
points  and  40  meV  at  nodal  points. 

However,  there  are  two  obvious  difficulties  with  this  scenario.  The  first  is  that 
the  spectral  peak  measured  by  ARPES  near  the  (ir,  0)  point  is  known  to  sharpen 
dramatically  when  one  goes  below  Tc,  a  phenomenon  interpreted  as  the  formation 
of  a  coherent  quasipartiele  in  the  superconducting  state.  This  sharpening  has 
normally  been  attributed  to  the  well-known  collapse  of  the  inelastic  scattering  rate 
below  Tc  due  to  the  opening  of  the  superconducting  gap,  but  it  is  hard  a  priori  to 
guess  why  something  similar  should  happen  in  the  presence  of  an  elastic  scattering 
rate  of  order  lOOmeV.  The  second  problem  is  that  recently  increased  momentum 
resolution  [34]  and  the  use  of  different  photon  energies  [102,  103,  104,  105]  has 
resolved  a  bilayer  splitting  which  has  its  maximum  effect  near  the  (7r,0)  point. 
Some  of  the  previously  observed  "elastic  broadening"  is  therefore  certainly  due  to 
this  as  well  as  to  pseudogap  effects  [106],  but  exactly  how  much  is  not  clear. 

On  the  other  hand,  one  can  certainly  not  ignore  the  out-of-plane  disorder.  The 
BSCCO  material  is  thought  to  be  doped  by  excess  oxygen  in  the  SrO  and  BiO 
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planes,  and  even  the  best  single  crystals  are  believed  to  contain  significant  amounts 
of  cation  switching  and  other  out-of-plane  defects  [107].  It  is  therefore  reasonable 
to  assume  that  quasiparticles  moving  in  the  Cu02  planes  of  this  material  must 
experience  a  smooth  potential  landscape  due  to  these  defects,  and  useful  to 
pursue  the  question  of  the  effect  of  this  type  of  scattering  in  the  superconducting 
state.  In  fact,  fits  [83]  to  Fourier  transformed-scanning  tunnelling  spectroscopy 
measurements  [69,  70]  on  similar  samples  to  those  used  in  the  ARPES  studies 
have  recently  been  shown  to  require  both  a  strong  (near-unitary  limit)  scattering 
component,  attributed  to  native  defects  in  the  Cu02  planes,  as  well  as  a  weaker, 
smooth  scattering  potential  component  attributed  to  defects  away  from  the  plane. 
Recently,  Markiewicz  has  also  attempted  to  relate  STM  and  ARPES  data  in  the 
superconducting  state  assuming  a  smooth  potential  [108]. 

5.2     Elastic  Forward  Scattering  in  Bi2Sr2CaCu2Q8 
5.2.1     Self-energy 

Normal  State.  Consider  a  model  system  including  elastic  scatters  of  finite 
range  K~l,  with  a  concentration  nh  Assuming  the  single  impurity  takes  the  form  3 
V{r)  =  V0e~Kr,  we  can  calculate  its  two  dimensional  Fourier  component: 

2irnV0 
Vkk'  =  ((k-k<)*  +  ^'  (5J) 

where  Vo  sets  the  strength  of  the  potential.  The  self-energy  in  Born  approximation 

is  then  defined  through  Eqn  5.6 

As  the  range  of  the  potential  k~'  increases,  the  scattering  of  a  quasiparticle 

from  k  to  k'  becomes  peaked  in  the  forward  direction.  As  shown  in  Fig.  5-2, 


3  There  is  a  consensus  that  quasiparticle  moving  in  the  Cu02  plane  experience  a 
smooth  potential  landscape  due  to  defects  whose  detailed  structure  is  unclear  yet. 
However,  a  concrete  impurity  model  is  needed  for  any  computations.  We  adopt  the 
exponential  form  rather  than  any  other  particular  structure  such  as  Yukawa  law 
for  the  sake  of  analytical  convenience.  Our  qualitative  picture  is  robust  against  any 
specific  form  of  forward  scatters. 
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Figure  5-2:  Geometry  for  the  forward  scattering  process  in  which  a  quasiparticle 
scatters  from  k  to  k'. 

when  k  is  close  to  k'  and  both  are  not  too  far  from  the  Fermi  surface,  we  may 
parameterize  them  as 


k 
k' 


kf  +  k_i_k_L 

k/.-  +  qn  +  k'Lk'± 


(5.8) 
(5.9) 


where  q  =  k  -  k'  is  the  momentum  transfer  and  qj  its  component  parallel  to  the 
Fermi  surface.  The  unit  vectors  k±  and  k'±  are  the  projections  of  k  and  k'  onto  the 
Fermi  surface,  respectively,  such  that,  e.g.,  ek,  =  vF(k')k'±.  The  imaginary  part  of 
the  retarded  self-energy  Eqn.  5.6  becomes 


£"(k,w) 


mVnitVo)2    f    dk  dfti 
■n.iK2Vg2     f  rf(/f| 


/ 


37TO,Vn2  1 


^^    [(*)Vr 


(5.10) 


Eqn.  5.10  shows  exphcitly  that  in  the  limit  of  small  re,  the  self-energy  becomes 
more  and  more  sharply  peaked  "on  the  mass  shell"  w  =  tk,  as  a  generic  feature 
of  long-range  potentials.  The  angular  dependence  of  the  self-energy  in  the  limit 
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of  K,  — »  0  can  be  verified  exactly  as  the  same  as  that  of  1/vp (k).  Furthermore, 
I  should  point  out  the  self-consistent  treatment  (by  requiring  E[C°]  — >  S[G])  is 
important  in  forward  scattering  case  since  it  eliminates  the  van  Hove  singularity  in 
the  spectrum,  which  may  account  for  the  absence  of  such  peculiarity  in  STM  and 
other  tunnelling  experiments  on  BSCCO,  while  in  the  case  of  point-like  scatters, 
the  self-consistency  (in  Born  approximation)  only  introduces  a  correction  of  the 
order  of  J1  /E2F  to  the  non-selfconsistent  result,  which  may  be  neglected. 

Superconducting  State.  In  the  supercouduting  state,  the  self-energy  is  approxi- 
mated similar  to  Eqn.  5.6  as: 

k' 

=      ^ZaTa  (5.11) 

a 

with  Nambu  components  Sa  defined  as  following; 

E0(k,w)    =    n,  ^|K(k,k')|2     2       "      A2  ■  (5-12) 

k' 

Ei(k,w)     =    n,]T|K(k.k')|2 


and 


Si(k,w)    =    -n,  ^|V'(k,k') 


uj2- 

-  p2   -  A2    ' 

ck'        "k' 
fk' 

2 

Ak' 

A2 


(5.14) 


k'  ~       CK 

In  the  limit  of  k  <<  1,  the  asymptotic  forms  of  self-energies  can  be  derived 
analytically,  which  read: 


Ej(k,w) 


r0(k,. 


1\Jix. 


-A2 


Y,  *»  x 


■l/y^ 


A; 


Wjp(1c) 


(5.15) 


Here  s„  =  \ui\,  -Aksgiiu;,  and  vsgnw^uj2  -  A2,  for  the  Nambu  components  a  = 
0, 1  and  3  respectively,  and  T0(k)  is  the  normal  state  scattering  rate.  E3  vanishes 
on  the  Fermi  surface  tk  =  0  in  this  limit.  Specifically,  when  the  momentum 
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k  is  close  to  the  Fermi  surface  and  the  energies  u  are  small,  such  that  |  (efc  ± 

y/ui2  -  A?k)/nvF(k)\  <  1,  the  self-energies  reduce  to  the  following: 


5j(k,w)   ~   -r0fo) -J1-.  (5.16) 


■*k 


isflM  -  r„(k,)  -^e^=  (5.i7) 

Sg(k,w)    M    0,  (5.18) 

but  are  strongly  suppressed  due  to  energy  conservation  when  \^ui2  -  Aj;  - 
ek|  becomes  greater  than  kvf,  as  one  may  observe  in  Fig.  5-34  .  The  detailed 
derivation  of  superconducting  self-energies  are  included  in  the  appendix.  It  is 
noteworthy  that  for  k's  away  from  the  node,  these  results  are  identical  to  those 
expected  from  an  s-wave  superconductor  [109]  (even  when  self-consistency  is 
included)  since  the  extreme  forward  scattering  eliminates  the  momentum  averaging 
over  the  Fermi  surface.  We  therefore  expect  to  recover  Anderson's  theorem,  the 
insensitivity  of  bulk  thermodynamic  properties  to  nonmagnetic  scattering.  The 
physics  here  is  that  the  for  k  away  from  the  nodal  direction,  forward  peaked 
scattering  doesn't  mix  the  order  parameter  of  different  signs,  and  hence  doesn't 
break  the  Cooper  pairs. 

In  Fig.  5  3,  I  also  present  the  numerical  results  for  the  self-energy  components 
at  nodal  and  antinodal  points  on  the  Fermi  surface.  The  bare  electron  band 


4  Roughly  speaking,  the  change  of  momentum  <5k  in  the  scattering  is  cut  off  by 
K.  Therefore  for  scattering  processes  close  to  the  Fermi  level,  the  corresponding 
variation  of  energy  is  KVpQs),  if  we  linearize  the  electron  dispersion.  On  one  hand, 
the  quasiparticle  energy  should  be  conserved  in  scattering  and  kinematics  only 
gives  a  finite  available  phase  space,  i.e.,  the  constant  energy  contour  as  discussed 
in  Chapter  3.  On  the  other  hand,  the  linearization  in  analytical  computation  may 
violate  the  exact  energy  conservation.  But  the  energy  difference  between  the  initial 
(ek)  and  final  electron  (y/w*  -  A£)  component  should  still  be  cut  off  as  well  due  to 
the  forward  nature  of  scattering. 
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Figure  5-3:  The  self-energy  terms  -Im  S0(k,w),  Im  Ei(k,a;),  and  -Im  2s(k,fa;)  in 
the  superconducting  state  at  T  =  0  for  k  =  kA  (top)  and  kw  (bottom),  for  k  =  5 
and  0.5  and  the  same  band  and  scattering  parameters  as  previously  used.  Here 
Ak  =  A0  (cos*  -  cos  ky)/2  with  A0  =  0.21. 
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dispersion  is  taken  as: 

tk  =  -2t(cos  kx  +  cos  ky)  -  it'  cos  kx  cos  ku  -  p.  (5.19) 

and  t/f  =  -0.35,  p/t  =  -1. 

The  quasiparticle  properties  near  the  Fermi  surface  iare  determined  by  some 
particular  combination  of  the  Nambu  self-energy  components,  by  examining  the 
denominator  of  the  full  Green's  function: 

^     _ek_Ak 

where  u)  =  lo  -  E0,  ek  =  £k  +  S3,  Ak  =  Ak  +  £,.  The  total  elastic  scattering  rate 
will  be: 

r,.,(k)  ~  -to  (s„(k.u;)  +  ^t  £i(k,«o)  ,  (5.21) 

provided  one  can  neglect  E3  (see  Fig.  5-3). 

In  the  forward  scattering  limit  where  Eqn.  5.16-5.18  hold,  the  effective  elastic 
scattering  rate  (  Eqn.  5.21)  becomes 


r,,(k.u;)~r1,(k^v/u;'      ^       ,      M>|Ak|.  (5.22) 

w\ 

For  k  along  the  nodal  direction,  the  elastic  broadening  in  the  superconducting  d- 

wave  state  is  equal  to  its  value  in  the  normal  state.  However,  for  k  in  the  antinodal 

region,  the  broadening  vanishes  as  u  -»  Ak  and  approaches  the  normal  state 

value  only  when  u  becomes  large  compared  with  Ak.  The  elastic  contribution 

Tei(kA,u>)  to  the  broadening  at  the  antinodal  kA  point  versus  u  is  shown  in 

Fig.  5-4.  Physically,  the  individual  contributions  to  the  normal  So(k,w)  and 

anomalous  E^kjUj)  self-energies  are  both  enhanced  by  the  density  of  states  factor 

(id  -  A2(k))~2  (Fig.  5-3).  However,  the  normal  contribution  describing  the 

scattering  out  of  state  k  into  k'  is  compensated  by  the  anomalous  contribution 

scattering  into  k  from  the  pair  condensate.  This  gives  rise  to  the  suppression  of  the 
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Figure  5-4:  Scattering  rate  I\,,(k,w)  vs.  u  for  k  =  k^left)  and  kjv(right)  in  the  su- 
perconducting state  at  T  =  0,  for  k  =  5  (top)  and  k  =  0.5  (bottom).  Here  r0(k  ,)  = 
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elastic  scattering  rate  seen  in  Fig.  5-4  relative  to  Fig.  5-3  as  u  approaches  -|AkJ 
from  below. 
5.2.2     Spectral  Function 

General  Spectrum.  In  the  forward  scattering  limit,  with  the  self-energy  given 
by  Eqns.  (5.16)-(5.18),  one  obtains  a  result  for  the  Green's  function  previously 
discussed  by  Markiewicz  [108], 

C(\c    A  ~  (uTo  +  &kn)z(k,u)  +  tkT3 

— lk-^)-— 7-5 A2,     ,, T5 B— •  5.23 

Here  z(k,u)  =  1  +  iToQ^sgaw/y/u2  -  A|.  The  electron  component  of  the  spectral 
function  is  then 

A(k,ui)    =    — ImGufow) 

71" 

1  T  u>zCk,w)  +  ek 

=  ^Im(^-A^(k,^-4-  (5'24) 

It  is  useful  to  consider  a  few  special  cases  of  (5.24)  more  closely.  In  particular,  on 
the  Fermi  surface  ck  =  0  one  has  the  simple  expression 


r0(k) i \oj\_ 

t    v/w2  -  A2  «/2  -  A2  +  r0(k)2 


*(kf  ,«„)  =     ^       r__^  ,.,2_J   '  ,2,  (5.25) 


while  near  the  gap  edge,  e.g.,  u  <  -|Ak|, 


A(k,W)«i    ^ M_ 

*  c2  +  r§(k)  v^^a2 


(5.26) 


At  the  nodal  point  kw,  where  the  gap  vanishes,  the  spectral  weight  is  given  by 
the  simple  Lorentzian  form 

</i         \  r0(kAi)/7r 

A(kN,u   =         '      n,    Y2-  (5-27) 

w2  +  r0(kw)2  v      ' 

and  at  low  temperatures  where  the  elastic  scattering  is  dominant  one  can  determine 
r0(kw).  However,  for  k  at  the  antiuodal  point  kA  such  that  -Ak4  -  6uj  <  ui  < 
~Ak„ 

1     Ak4  1 
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where  the  "width"  Su  depends  upon  the  ratio  AkA/Fa(kA).  If  Akl  is  large  com- 
pared with  Va(k.A),  the  width  Su  Si  r0(kA)2/(2Ak/1).  If  one  integrates  A(kA,u>) 
from  -Ak„  -  6u  to  -Ak/(  to  define  a  "peak  intensity" 
f- At 


/(k,,) 


f-a*A  1 

:  /  A(kA,uj)du)~  -,  (5.29) 


which  is  independent  of  Ak/,/r0(k,i).  However,  when  the  system  is  sufficiently 
dirty  such  that  Ak/1  <  Tofk^),  the  falloff  of  A(kA,ui)  as  u  decreases  below  -Ak, 
varies  as  (J2  -  Aj^)-1/2.  In  this  case,  the  scale  is  set  by  Ak/l  and  if  one  takes 
Su  =  Akj,  the  peak  intensity  varies  as  I(kA)  ~  Ak/l/ro(k^),  This  is  quite  different 
from  the  usual  BCS  quasiparticle  result  which  is  proportional  to  the  quasiparticle 
renormalization  factor  z(kA)  times  a  coherence  factor  which  is  1/2  on  the  fermi 
surface.  It  should  be  possible  to  test  the  foward  scattering  scenerio  by  comparing 
the  variation  of  I{kA)  with  Akj/ro(k,4). 

Nodal  and  Antinodal  Spectra.  To  compute  a  realistic  spectral  function  and 
compare  it  with  experiments,  one  should  include  a  small  concentration  of  unitary 
scatters  as  well  as  inelastic  electron-electron  collisions.  The  former  is  verified 
by  STM  experiments  and  has  been  intensively  discussed  in  Chapter  2-3.  With  a 
concentration  roughly  ~  0.2%  decided  from  the  number  of  zero  bias  resonances  in 
LDOS  map,  it  gives  a  broadening  in  the  self-consistent  T-matrix  approximation: 

(5.30) 


—■»       £kc(k,^)'u- 

This  gives  a  scattering  rate  of  T,,  ~  10~3(  in  the  normal  state,  leading  to  an 


impurity  bandwidth  7„  ~  VruA0  ~  W~2t  (of  order  1  to  2  meV).  As  far  as 
the  inelastic  self-energies  are  concerned,  we  use  a  numerical  interpolation  of 
the  u—  and  T-  dependence  obtained  from  the  spin-fluctuations  performed  by 
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[110,  111,  112].5    Combining  the  contributions  from  these  three  channels,  the  total 
scattering  rate  is  then: 

£,o,=S,(,/ +£.,!,„  +Sinel-  (5.31) 

The  spectral  functions  at  antinodal  point  A  and  nodal  point  N  with  the 
scattering  rate  defined  in  Eqn.  5.31  and  the  Fermi  function  cutoff  }(ui/T)  are 
plotted  in  Fig.  5-5.  In  the  antinodal  direction,  a  square-root-like  behavior  develops 
when  u  approaches  -Akj  from  below.  This  should  be  contrasted  with  the  broad 
Lorentizian  peak  in  the  normal  state  and  is  an  intrinsic  consequence  of  elastic 
forward  scattering.  It  is  also  noteworthy  to  point  out  that  this  anomaly  at  Ak  is 
stable  against  the  many-body  interactions  since  at  lower  temperature,  when  u)  is 
equal  to  the  real  part  of  the  gap  at  the  gap  edge,  the  inelastic  broadening  vanishes 
as  T3  and  the  elastic  components  dominate. 

Along  the  nodal  direction,  since  the  scattering  always  mixes  the  initial  and 
final  states  with  order  parameters  of  different  signs,  a  Lorentizian  centered  at  Fermi 
level  is  observed  for  zero  temperature.  This  peak  will  be  further  cut  off  by  the 
Fermi  distribution  and  broadened  by  inelastic  effects  when  temperature  increases. 
In  the  limit  of  k  =  0,  we  recover  the  normal  state  spectral  function,  even  at  zero 
temperature 

Quasiparticle  Dispersion  Near  the  Antinodal  Point.  In  a  clean  superconductor, 
there  is  a  peak  in  the  spectral  function  A(k,w)  at  the  quasiparticle  pole  u  =  £k  = 
\Ak  +  Ak-  ln  particular,  as  the  momentum  moves  along  the  cut  from  (kA,  jr)  to 
(0, 7r),  one  expects  to  see  a  dispersion  of  this  peak  to  higher  energies. 


5  The  spin-fluctuation  self-energies  behave  like  follows:  at  the  nodal  point,  at  low 
temperatures,  the  scattering  rate  initially  increased  as  the  third  power  of  w  or  T 
depending  upon  which  is  larger.  At  other  k  points  on  the  Fermi  surface,  the  scat- 
tering rate  varies  approximately  as  the  third  power  of  this  energy  measured  relative 
to  Ak.  The  reduction  of  the  inelastic  scattering  rate  at  low  excitation  energies  re- 
flects the  suppression  of  the  low  energy  spin-fluctuations  due  to  the  opening  of  the 
d-wave  gap. 
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Figure  5-5:  Finite  temperature  spectral  function  at  the  antinodal  point  A  and 
nodal  point  N  on  the  Fermi  surface  multiplied  by  the  fermi  function.  Results  for 
k .  =  2  and  0.5  with  T0(kA)  =  0.2t  are  shown. 
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Figure  5-6:  A(k,w)  vs.  u)  for  k      =       2  and  0.5.  Results  are  given  for  the  k  points 
at  (tt,  0),  {it,  0.057r),  (it,  0.l7r),  (it,  0.157r).  The  disorder  levels  correspond  to 
r0{kA)/A0     =     1,  0.5,  and  0.025.  Note  the  spectra  for  different  k  points  have  been 
offset  for  clarity. 


However,  if  the  forward  elastic  scattering  strength  ra(kA)  >  AkA,  then  the 
peak  in  .4(k,w)  remains  at.  -At  rather  than  dispersing.  Fig.  5-6  shows  plots  of 
A(k,uj)  for  different  values  of  k  between  the  M  and  A  points  for  n  =  2  and  0.5,  and 
several  values  of  the  scattering  rate  ro(k.4). 

As  samples  improve,  there  is  a  natural  tendency  in  this  model  for  the  spectrum 
for  k  not  too  far  from  the  antinode  to  cross  over  from  one  characterized  by  a 
nondispersive  peak  at  Ai<  in  the  dirty  limit  where  T0  ~  Ao  to  one  characterized 
by  a  dispersive  quasiparticle  peak  at  Ey,  when  T0  is  small  compared  to  A0.  This 
crossover  is  due  to  the  way  in  which  the  forward  elastic  scattering  rate  for  a  d-wave 
superconductor  is  reduced  as  the  gap  edge  is  approached  and  is  analogous  to  the 
same  effect  discussed  analytically  (in  Appendix  B)  for  an  s-wave  superconductor. 
In  a  system  with  T0  >  A0,  no  quasiparticle  peak  is  observed,  but  a  sharp  feature 
does  appear  at  -Ak,  representing  simply  the  spectral  weight  in  the  overdamped 
quasiparticle  peak  piling  up  at  the  gap  edge  as  in  the  s-wave  case.  Only  when  To 
becomes  small  compared  to  A0  does  one  see  a  true  quasiparticle  peak  dispersing 
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Figure  5-7:  Comparison  between  recent  ARPES  data  (left  panel)  and  calculated 
A(k,ui)  (right  panel). 


as  —Ek.  In  the  most  strongly  forward  scattering  case,  re  =  0.5,  one  can  see  that, 
depending  on  the  strength  of  the  scattering  rate,  one  can  have  simultaneously  a 
broadened  dispersing  feature  as  well  as  a  gap  edge  feature. 

In  Fig.  5-7,  the  measured  ARPES  spectra  along  a  special  cut  (  (0,  0.64jt)  — > 
(kf,  0.647r)  )  (data  from  private  communication  with  Stanford  ARPES  experi- 
mental group  (2004))  and  the  theoretically  calculated  spectra  adopting  elastic 
impurities  with  8%  and  re  =  0.5  are  compared.  While  the  variation  of  the  calculated 
order  parameter  is  smaller  than  the  exact  measurement,  other  details,  such  as  the 
dispersion  of  the  hump  feature  and  the  gap  peak  intensity,  are  in  good  agreement 
with  each  other. 

5.3     Discussion  and  Conclusion 

Elastic  forward  scattering  of  quasiparticles  due  to  spatially  smooth  potentials 
is  prevailing  in  BSCCO.  We  have  discussed  how  this  scattering  can  fit  with  the 
BCS  framework  to  invalidate  out  naive  intuition  on  the  insensitivity  of  elastic 
scattering  self-energy  to  entering  the  SC  state  and  hence  argued  that  elastic 
forward  scattering  does  not  necessarily  prevent  the  formation  of  quasiparticles 
across  the  normal-superconducting  transition. 
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While  the  bulk  and  thermodynamics  properties  such  like  Tc  is  [113]  and  the 
transport  properties  such  as  conductivity  are  not  significantly  modified  by  extended 
impurities  because  of  their  forward  nature(there  is  not  significant  momentum 
transfer  in  scattering),  the  spectroscopic  properties  are  modified  significantly  by 
the  elastic  forward  scattering:  the  near-cancellation  of  the  two  Nambu  components 
of  the  self-energy  near  the  gap  edge  in  the  forward  scattering  limit  leads  to  a 
dramatically  reduced  elastic  scattering  rate  in  the  superconducting  state,  which 
sharpens  the  spectral  features  of  quasiparticles  which  are  not  too  close  to  the 
nodes.  Besides  the  collapse  of  self-energies,  a  square-root-like  singularity  local  is 
also  found  to  develop  at  the  local  gap  edge  u  =  Ak. 

The  local  gap  anomaly  is  of  particular  interest  since  it  is  inherent  from  forward 
scattering.  This  singular  structure  remains  when  k  is  below  the  Fermi  surface 
and  is  speculated  to  be  robust  against  many-body  effect  because  the  latter,  in 
general,  is  suppressed  when  temperature  decreases.  To  date,  ARPES  experiments 
generally  mapped  out  the  momentum  anisotropy  of  Ak  along  the  Fermi  surface, 
with  assumed  information  of  bare  electron  dispersion  c*.  However,  this  gap  edge 
anomaly  seems  to  propose  an  alternative  method  to  determine  the  gap  function  for 
k  values  away  from  Fermi  surface  by  taking  use  of  forward  scatterers.  Successful 
realization  of  this  technique  may  complete  the  momentum  dependence  of  Ak  in 
the  whole  first  Brillouin  zone,  which  definitely  contributes  to  the  understanding  of 
cuprates.  Furthermore,  another  possible  advantage  of  this  method  is  that  it  doesn't 
necessitate  any  knowledge  of  ej,  in  advance  at  all! 

Moreover,  in  samples  of  less  disorder,  we  found  a  dispersive  peak  associated 
with  the  bare  quasiparticle  dispersion.  It  is  tempting  to  speculate  that  this  phe- 
nomenon is  related  to  the  peak-dip-hump  features  observed  generically  below  Tc  in 
cuprate  ARPES  experiments.  When  the  system  becomes  dirtier,  the  quasiparticle 
peak  will  be  broadened  from  <5-like  function  into  Lorenzian  form.  However,  far  away 
enough  from  the  nodal  direction,  since  we  physically  recover  Anderson's  theorem 
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and  Ak  remains  unchanged,  the  weight  of  the  Lorenzian  tail  cannot  leak  into  the 
subgap  regime,  namely,  \uj\  <  |Ak|  and  will  keep  piling  up  at  the  leading  gap 
edge  instead.  We  then  expect  this  anomaly  is  qualitatively  more  pronounced  in 
dirties  samples.  This  could  possibly  explain  the  excessive  peak  height  observed  by 
STM  experiment  at  the  coherent  gap  edge,  which  is  even  much  higher  than  the 
theoretically  calculated  peak  height  for  pure  systems. 

The  last  concern  of  this  chapter  is  the  lifetime  of  nodal  quasiparticles.  By 
extracting  the  half-width  Ak  of  the  peaks  in  MDC,  we  obtained  a  scattering 
rate  which  is  linear  in  temperature  above  Tc  but  varies  nonhnearly  below  Tc  and 
saturates  to  a  value  determined  by  the  small-angle  scattering.  However,  it  was 
reported  by  Valla  et  al.  [16]  that  the  width  of  MDC  has  a  linear  dependence  on 
temperature  across  the  normal-superconducting  transition,  which  suggested  the 
irrelevance  of  the  nodal  quasiparticles  in  the  superconductivity  and  a  quantum 
critical  behavior.  This  discrepancy  is  rather  interesting  and  definitely  worthy  of 
further  investigation. 


CHAPTER  6 
CONCLUSION  AND  REMAINING  QUESTIONS 

I  have  studied  impurity  induced  quasiparticle  interference  and  the  related 
electronic  structure  in  fully  disordered  BSCCO  samples.  I  adopted  two  types  of 
nonmagnetic  impurities:  short-ranged,  unitary  scatters  and  spatially  extended, 
weak  scatters,  in  conjunction  with  the  d-wave  version  of  the  mean  field  BCS 
theory  (which  is  believed  to  be  appropriate  to  optimally  to  overdoped  cuprates),  to 
calculate  the  spectroscopic  properties  such  as  local  density  of  states,  local  gap  map 
and  spectral  function. 

I  discussed  how  the  theoretical  predictions  on  Zn  and  Ni  impurities  mesh  with 
the  details  of  the  atomic  scale  resonances  observed  in  the  STM  image  of  BSCCO. 
I  addressed  on  the  dispute  regarding  the  characteristics  of  Zn,  which  is  not  fully 
pinned  down  yet. 

There  are  divergent  theories  concerning  how  quantum  interference  between 
many  impurities  modifies  the  low  energy  quasiparticle  states  and  the  Fourier 
transformed  DOS.  I  joined  into  this  discussion  by  solving  the  two  impurity  problem 
exactly.  I  concluded  that  the  divergent  p(uj)  at  the  Fermi  level  for  half-filled  nearest 
neighbor  band  is  the  consequence  of  the  special  r2  symmetry  of  the  nearest  tight- 
binding  model  seen  in  the  analogous  the  two  impurity  interference  effect.  I  further 
proposed  a  real  space  perspective  regarding  the  homogeneous  broadening  by  finite 
STM  resolution  and  explained  why  the  impurities  appear  nearly  isolated  in  the 
STM  image  while  they  are  expected  theoretically  to  interfere  with  each  other  even 
separated  by  long  distance. 

I  further  examined  the  "octet"  single  impurity  model  as  a  way  understanding 
the  Fourier  transformed  DOS  in  BSCCO.  I  argued  that  while  its  cm-shell  approx- 
imation captures  the  fundamental  scattering  processes,  it  fails  to  explain  other 
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details  such  as  the  peak  widths  and  peak  height  for  certain  scattering  wavevectors. 
I  showed  why  it  is  necessary  to  include  both  of  the  two  kinds  of  disorder  above  to 
reconcile  these  discrepancies.  Additionally,  I  discussed  the  implications  of  extended 
impurities  to  the  nanoscale  gap  inhomogeneities  in  STM  measurements. 

The  extended  weak  potentials  have  a  remarkable  influence  on  the  spectral 
function.  From  my  computation,  I  realized  that  the  quasiparticle  scattering  rate 
actually  collapses  in  the  superconducting  state  due  to  the  near  cancellation  between 
the  normal  and  anomalous  scattering  processes,  yielding  well-defined  antinodal 
quasipartieles.  This  interesting  finding  invalidates  the  argument  that  elastic 
scattering  is  intuitively  robust  against  the  opening  of  a  gap  upon  lowering  the 
temperature. 

In  my  dissertation,  I  have  established,  within  the  frame  of  quasiparticle 
scattering  due  to  impurities  of  different  spatial  nature,  a  reasonably  coherent 
picture  subsuming  recent  observations  in  STM  and  ARPES.  However,  despite  this 
modest  success,  there  are  still  remaining  discrepancies  between  STM  and  ARPES 
experiments  from  this  point  of  view,  such  as,  1)  why  we  have  nanoscale  electronic 
disorders  at  the  gap  edge  but  quasipartieles  are  well  defined  in  momentum  space? 
2)  why  the  ARPES  spectra  for  the  nodal  quasiparticle  display  a  kink  while  STM 
does  not  reflect  any  collective  mode  features.  Full  comprehension  of  these  puzzles 
may  require  physics  beyond  the  quaisparticle  scattering  scenario.  However,  I  would 
like  to  put  a  period  here,  as  I  realized  that  I've  been  too  talkative. 


APPENDIX  A 
GREEN'S  FUNCTION  FOR  THE  SUPERCONDUCTING  STATE 

In  this  section,  we  try  to  evaluate  the  real  space  component  of  the  Green's 

function  for  a  pure  nested  band,  i.e.,  ek  =  -2l(coskx  +  cosfc,).  If  we  define, 

r  =  (m,  n),  where  m,  n  are  in  unit  of  lattice  constant  a,  the  real  space  Green's 

function  is  defined  through  the  Fourier  transform  as: 

G°(r,w)    =    ^eik'rG°(k,w), 
k 

=    2_,  cos(kxm)  cos(fc„n) 

k 

wfo  +  (e*,  -  n)f3  +  Akfi 


"2~El 


(A.l) 


Start  from  Eq.  (A.l),  and  we  express 


cos(fcxm)    =    2"*  '  cos™  kT 

+  2j^    l)     j[(m-2j)\    ^cos/^ 

where  [. . .]  refers  to  the  integer  part  of  the  argument.  We  focus  on  the  half-filled 
case  ji  =  0  and  write  Eq.  (A.l)  as  the  sum  of  terms  of  the  form 

f-lfi    =     >™s'(fcr)cos''(A:!, — 

k  "2  -  H 

where  p  =  m,  m  -  2, . . .  and  q  =  n,  n  —  2 We  proceed  by  linearizing  the 

dispersion  near  the  node  at  (7r/2, 7r/2)  and  making  the  coordinate  transformation 
£2  =  4  +  Ak,  tan<?  =  Ak/ck. 

_o_   f2"  dO  f    smfl      cos eV  /sine      cosO 


2J>+«  yD     2tt  V      A0  2/    J    \  A0  2( 

.  [A i^^+^dE—"  +  E^os ^  +  sm "fl ' 
Jo 


w2-£2 
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The  prefactor  is  a  =  N/(2irvFv&)  where  N  =  A  is  the  number  of  nodes,  vF  is 
the  Fermi  velocity  and  uA  is  the  anomalous  quasiparticle  velocity  |VkAk|,  and  the 
cutoff  A  is  of  order  A0.  The  integrals  over  E  and  9  are  easily  done  and 

where  P^  are  constants  given  by  the  angular  integrations,  and 


FaM      =        f 

Jo 


fjdr; 


E" 


i-y  - 


The  constants  P3m  vanish  for  j  =  1,3  when  p  +  q  =  even  and  vanish  for  j  =  0  when 
p  +  q  =  odd.  The  first  few  nonzero  elements  are 


p°    - 

1 

P1    - 

Mo     — 

01~2A0 

P3     - 
Mo     — 

'm~     At 

P°     - 

Ml      — 

11          2A§  +  8t2 

Only  even  moments  of  Fa(uj)  are  needed: 

Since  we  are  interested  in  the  leading  order  behavior  of  (7(R,  w)  we  note  that  for 

small  ui  , 


1       A2 

2h'^ 

A2^ 

2n 


For  R  =  (2m,  2n),  the  leading  order  contribution  to  G(R,uj)  comes  from  the 
single  term  in  the  expansion  containing  gm.  To  second  order  in  ui: 


(?(R,w)  =  -(-l)"""^ln^f0  +  wCo(R)f„, 


(A.2) 
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where  C0(R)  is  real,  and  is  the  sum  of  several  terms.  The  largest  term  contributing 

to  C(R)  is  of  order 

q|v| 


n)  \A0) 


2(m+n) 


16(m  + 
from  which  we  estimate  a  range  of  validity 

(A/A„)'CH"> 

M  <  Ae       "=<"•+">     . 

For  other  R,  there  is  no  single  dominant  term  in  the  expansion  for  the  Green's 
function,  and  the  leading  order  behavior  comes  from  the  sum  over  a  large  number 
of  real  nondivergent  terms.  For  our  purposes,  it  is  sufficient  to  note  that  when 
R  =  (2m  +  1,  In  +  1),  the  sums  take  the  form 


G(R,w)  =  wC0(R)f0, 

and  when  R  =  (2m.  +  1, 2n)  or  (2m,  2n  +  1) 

e?(R,<j)=Ci(R)fi  +  Cs(R)f&, 

where  C0(R),  Ci(R),  and  C3(R)  are  real  constants. 


(A.3) 


(A.4) 


APPENDIX  B 
EVALUATION  OF  THE  SELF-ENERGY  IN  BORN  APPROXIMATION 

In  this  section,  I  give  out  the  evaluation  of  self-energy  in  Born  approximation 

for  several  different  types  of  impurities. 

B.l     "Yukawa  I"  Potential 

Assume  that  V(r)  =  Qexp(-Kx)/r,  and  its  Fourier  transform  in  2D  takes  the 


form  Vq  =  2ttQ/v/(/2  +  k2.  According  to  Fig.  5-2,  we  obtain: 

k    =    kr  +  kj;±  (B.l) 

k'    =    kp  +  qg  +  fc^  (B.2) 

for  momentum  close  to  the  Fermi  surface,  we  can  linearize  the  dispersion  and  we 
have  tk'  =  VF(k')k'±. 
B.l.l     Normal  State 

The  self  energy  in  Born  approximation  can  be  then  evaluated: 
-u^ttQ)2 


'  (2 


2irQY      f      ,    f      ,  1 

^*  J  *l  J  ^(k-kV  +  ^  %  "  *) 

=    -nmCffdE.     *     , 5 ! , 

M   _Inj(l_   r      i 

\vF(h)\J       ■  [q||'+fcx(fcx__^y)]2  +  AcS 
™,Q2      /•     ,  I 

7r2njQ2 

vWMfc|iH?+ARV'  (B'3) 

where  in  the  middle  steps  the  difference  between  vF(k')  and  vF{k)  is  neglected 
since  in  the  limit  of  small  k  the  two  approach  one  another,  and  we  need  to  retain 
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only  the  most  singular  dependence  on  k-k'.  In  addition  it  is  reasonable  to  set  that 
dk\\  =  d</||,  since  k  is  a  constant  in  the  integration. 
Special  Cases:.  •  k  =  k>*  (on  Fermi  surface) 

S"(kF,c)  =  -^=T  (B.4) 

•  k  =  kp,  u  =  0 

S"(kF,w  =  0)  =  -i^-.  (B.5) 

Note  this  diverges  in  forward  scattering  limit  as  k  — >  0.  If  we  normalize  to  EJ,  this 

gives 

S"(V,u,  =  0)       |MM 

s"(k4,^  =  o)    Mk)r  ,u-0' 

•  u)  =  0,  k  7^  k/? 

S"(kF,  w) 7r2"";Q'  (B.7) 

Note  in  forward  scattering  limit  k  — »  0,  self-energy  is  a  very  singular  function  of 

either  k±  or  id. 

B.1.2     Superconducting  State 

In  the  main  Nambu  channels,  i.e.,  tauo,  taui,tau3,  the  imaginary  parts  of  the 
Green's  function  are  the  summation  of  two  5  functions: 


ImG0(k',u;)    =    Im      _     2   =  --  [S{u  -  Ek.)  +  5(ui  +  Ek,)}  (B.8) 

Imdfk'^)    =    Im   2Akr2   =-l^[6(u-Ek,)  +  &{u,  +  Ek,)}  (B.9) 

lmG,(k',L0)    =    Im       'k        =  ~(±  [S{uj  -  Ek.)  +  5(lu  +  Ek,)j ,  (B.10) 

W     —  Pj^i  Z  U) 
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where  Ek,  ~  JvF(k'u)2k'±2  +  A£,.  The  S  functions  in  the  above  equations  can  be 
further  evaluated  as  follows: 


S(u,±Bv)    ta    J^ 


oi-k- 


ey. 


-°^S?i 


(B.ll) 


^k'L2  +  Ak, 


upi 


S\k±- a     ...  ,,.a       I  «(Tw)         (B.12) 


"F(k') 


13) 


The  sum  of  2  delta-functions  is  therefore 
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^  %  +  a£k.) 


^rE/k-^?V(-) 
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"f(k') 


.V^3^ 


rp  Component  of  the  Self-energy. 


Eo(k,i<;) 


n,OV 

■K 

l"F(k)| 

2 

n.QV 

7T 

<-F(k) 

2 

»//-'k' 

jr 

|t',.(k)|  2 


jdk[Jdk'L 


E^Mfci-a^£ 


qii  +  kxkj.  -  k'±k'± 


(V^-Al,)-' 


m  +  [k 


«,^U<, 


(v^3^)-1 


_±  «|  +  [k±-cui 


.(B.14) 


The  integral  as  it  stands  is  not  trivial  due  to  the  dependence  of  Ak-  on  q^.  Since 
k'  is  close  to  k',  gy  is  small  so  we  may  write  Ak<  ~  Ak  +  v&q%.  The  most  singular 
factor  in  the  limit  k  ->  0  is  the  {k2  +  a'2)-1  factor,  so  we  are  justified  in  neglecting 
the  0||  dependence  in  the  square  root.  The  self-energy  then  takes  the  form: 
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£o'(k,ui) 


h,Q2tt\uj\ 


j  din\  Y, 


l  +  fa-ayffiC)   +* 


7r2n,Q2        \u\ 


2       V^A* 

x     Z     . 

—±  V(fcxt>H*||)  -  «vA>2-A2)2  +  «2t»f  (fc||)a 

Obviously  the  result  is  valid  only  for  u  >  Ak;  otherwise  the  Im  part  of  E  is 
On  the  Fermi  surface,  with  k  =  k,, ,  this  reduces  to 


(B.15) 


(B.16) 


E"(kF,w)~-7rV<22 


M 


Note  this  is  NOT  the  form  we  had  earlier  in  Chapter  5. 


(B.17) 


E'0'(k,.,W)~-r(k)M/v/^-A2,  (B.18) 

unless  u>2  -  Ak  <<  K2)>F(jfc||)2. 

t3  Component  of  The  Self-energy.  For  the  t3  component,  proceeding  similarly, 


we  get 


E3(k,w) 


^   ^  '  lr.   J-t.t.     _   i-'   J.'     1 


"  Mk)| 

n^sgnwjr 


(q||  +  fcj_A-x  —  *l*:l)    +k2 


qn+  (fci.-W1^''  )^A 


jt^Q2  ^ 
—  sgnw^- 


Here  the  presence  of  the  two  contributions  for  the  two  different  signs  of  k±  is 
crucial,  as  we  see  that  on  the  Fermi  surface,  k±  =  0,  E3  vanishes. 


(B.19) 
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Ti  component  of  the  self-energy, 

Ei(k,w)    ~     sgnw — 

x  y;  (b.20) 

<*=±  V(*x»f(*||)  -  aJuZ-AZf  +  «V(fc||)2 
Normal  State  Limit.  It  is  instructive  to  take  the  limit  Ak  — >  0  of  the  above 
expressions  and  attempt  to  recover  Eqn.  B.3.  E0  and  E3  become 

2  Aj<?2  y^  1 

2       ^i  ^(k±vr(klt)  -  a\u>\)*  +  K2vF(k^ 
y    _       n'2njQ2sg,nui  s~-y  a 

2    ~ hi  ^WMh)  -  «M)2  +  «2Mfc|i)2 '        (R22) 

Now  the  normal  state  self  energy  is  the  Nambu  11  component,  i.e.,  E0  +  E3,  which 
goes  to 


Note  odd  frequency  dependence  of  E3  is  necessary  to  recover  the  correct  result. 
B.2     "Exponential"  Potential 
Consider 

V(r)/Va  =  e-Kr  (B.24) 

which  gives  in  momentum  space  (neglecting  the  periodicity) 

Vic*    =    V0  f  rdrd<t>  e-^e^"** 

=    2-nV0  I     rdre-*'V0Mqr) 
Jo 


2nnV0  2ttkV0 


((/2  +  „2)3/2  ((k_k/)2  +  K2)3/2- 


(B.25) 
(B.26) 


B.2.1     Normal  State 

Self  energy  in  Born  approximation: 

S"(k,o;)    =    -  3,'iKy 


16|,.F(fc|,)| 


1  Special  ease  k  =  kp ,u  =  0: 

£"(kF,«j  =  0) 


(*x-^)a  +  «' 


37r2njK02    1 


.5/2 
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(B.27) 


(B.28) 


Clearly  the  above  form  is  more  singular  than  Yukawa  I  potential,  but  still  varies  as 
|tijr(k)|_1  around  the  Fermi  surface  for  any  nonzero  k. 
B.2.2     Superconducting  State 

If  we  are  allowed  to  perform  analogous  operations  as  we  did  hi  the  last  step 
and  neglect  the  q  dependence  of  the  gap,  we  find 


Ei'(k,u)    a 
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Trni^Vo2!^! 
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3ir2niK2V'a2  r— > 


«i  +  (*x-«V^r)  +  *2 
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(B.30) 


(B.31) 


B.3     Delta-function  Type  Potentials 
What  type  of  potential  gives  a  self  energy  as  in  Chapter  5,  e.g., 


[]"(k,v«>  =  o)  ~  r0(k)H/^2-Aj 


(B.32) 
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It's  not  a  pure  2D  delta  function.  If  we  simply  assume  \Vkk,\'2  =  Vo^k.k',  the  self 
energy  in  the  superconducting  state  is 

ES(k,w)    =    niV0^Ti-G"(k,u)  (B.33) 

=    "^riV-^-Ag  (B'34' 

=    -niV02^[S(ui-Ek)  +  6(io  +  Ek)},  (B.35) 

which  is  much  more  singular  on  the  Fermi  surface  than  Eqn.  B.32.  To  recover  the 
smoother  result  we  actually  need  to  assume 

|VW|2  =  V^Vj  (B.36) 

and  then  we  get  (see  Eqn.  B.13): 

4rr|«F(k)|  v^^A^'  (B-37) 

It  would  be  nice  to  ask  what  potential  would  result  in  real  space  from  a  ID  delta 
function  like  this.  One  way  to  regularize  the  delta  function  would  be  to  write 


vw  =  V(*n-$  +  "a  (R38) 

but  of  course  we  see  immediately  that  this  is  not  a  function  of  |k-k'|  and  therefore 
can't  correspond  mathematically  to  a  function  of  r  if  we  reverse  Fourier  transform. 
The  problem  is  that  to  get  the  special  form  of  the  self-energy  near  the  Fermi 
surface,  we  oidy  need  the  k-space  form  of  the  matrix  element  near  the  Fermi 
surface,  but  this  does  not  suffice  to  determine  the  form  in  real  space  everywhere. 
We  would  need  to  guess  a  real  space  form  which  gives  this  result  in  k  space  for 
a  given  band,  which  seems  difficult.  Therefore  the  only  connection  one  can  make 


109 


with  the  simple  form  in  Chapter  5  is  an  approximate  one.  Not  in  the  full  forward 
scattering  limit  k  — >  0  but  in  an  intermediate  regime  given  by 


(*l 


V^^l 


<£««  1 


(B.39) 


do  we  actually  get  exactly  the  anticipated  result.  It  is  interesting  to  note  that 
this  works  actually  for  any  form  of  the  scattering  potential.  The  only  thing  that 
changes  is  the  r0(k)  describing  the  normal  state  scattering. 
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